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Abstract. The experimental revolution in ultracold atomic gas physics over 
the past decades have brought tremendous amounts of new insight to the world of 
degenerate quantum systems. Here we compare and constrast the developments 
of cold atomic gases with the physics of nuclei since many concepts, techniques, 
and nomenclatures are common to both fields. However, nuclei are finite systems 
with interactions that are typically much more complicated than those of ultracold 
atomic gases. The simularitics and differences must therefore be carefully 
addressed for a meaningful comparison and to facilitate fruitful crossdisciplinary 
activity. We first consider condensates of bosonic and paired systems of fermionic 
particles with the mean-field description but take great care to point out potential 
problems in the limit of small particle numbers. Along the way we review some of 
the basic results of BEC and BCS theory, as well as the BCS-BEC crossover 
and the Fermi gas in the unitarity limit, all within the context of ultracold 
atoms. Subsequently, we consider the specific example of an atomic Fermi gas 
from a nuclear physics perspective, comparing degrees of freedom, interactions, 
and relevant length and energy scales of cold atoms and nuclei. Next we address 
some attempts in nuclear physics to transfer the concepts of condensates in nuclei 
that can in principle be built from bosonic alpha-particle constituents. We also 
consider Efimov physics, a prime example of nuclear physics transfered to cold 
atoms, and consider which systems are more likely to show interesting bound 
state spectra. Finally, we address some recent studies of the BCS-BEC crossover 
in light nuclei and compare them to the concepts used in ultracold atomic gases. 
While many-body concepts such as BEC or BCS states are applicable in both 
subficlds, we find that the interactions and finite particle numbers in nuclei can 
obscure the clear meaning they have in cold atoms. On the other hand, universal 
results from atomic physics should have impact in certain limits of the nuclear 
domain. In particular, with advances in the trapping of few-body atomic systems 
we expect a more direct exchange of ideas and results. 
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1. Introduction 

The last two decades have witnessed dramatic developments in terms of realizing 
Bose-Einstein condensates and degenerate Fermi gases which has been driven by 
the experimental advanced in trapping and cooling dilute atomic gases (Pethick & 
Smith 2008, Pitaevskii & Stringari 2003, Leggett 2006, Bloch et al. 2008, Giorgini 
et al. 2008). Typically, these system contain thousands or even millions of particles, 
but very recently samples with very small particle numbers have been realized 
(Serwane et al. 2011, Ziirn et al. 2012, Bakr et al. 2010, Sherson et al. 2010, Weitenberg 
et al. 2011). This regime is particularly interesting since the physics that can be 
learned here might aid our understanding of mesoscopic and microscopic systems in 
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other subficlds. One such venue is nuclear physics, where particle numbers are finite. 
This raises an interesting question of how concepts from the large system limit, such 
as condensation or superfludity, transfer to small systems. In the current presentation 
we will address similarities and differences between the physics of nuclei and ultracold 
atomic gases. 

There are a number of immediate contrasts that must be highlighted in a 
comparison of nuclear and cold atoms. First, the notion of scattering length and 
effective range expansion (Landau & Lifshitz 1981) is very efficient for cold gases since 
they are valid for very low energy and the typical collision energies approach the zero- 
energy limit. In more traditional atomic collision physics and in charged Fermi liquids 
in metals the long-range Coulomb interaction complicates this simple description of the 
two-body interaction. Secondly, nuclei are self-bound system implying that repulsive 
and attractive parts of the nuclear interaction must balance. In the atomic gas with 
external confinement there is no such requirement and the system can be studied in a 
wider range of regimes in principle. The relatively small particle numbers and the self- 
bound structure can have important consequences to the conceptual transfer between 
the subfields and will be a prime concern in this presentation. Furthermore, nuclear 
methods deal with structures resembling those of present interest in cold atomic gas 
investigations, i.e. (i) two inherently different but similar components, spin-up and 
down of neutrons and protons, as in mixed condensates, (ii) "unbound" boson pairs 
of nucleons in the supcrfluid BCS low-energy regime, (iii) boson expansions in terms 
of pairs of fermionic nucleon as molecular condensates, (iv) bosonic a-particle cluster 
structure as in one-component Bose condensates, (v) a-clusters and nucleons in mixed 
boson-fermion structures. 

In broad terms, nuclear and atomic-molecular physics have a lot of concepts, 
phenomena, and techniques in common. The reasons are to be found in the short-range 
character of the different interactions and the finite spatial extension of systems with 
a finite number of particles. Specifically characteristic units of length L (interparticlc 
spacing) and mass m can be used to construct a unit of energy E = h 2 / (jtiLq) . Using 
Lq and Eq as length and energy units provide dimensionless quantities which remove 
superficial large differences, emphasize similarities and allow focus on the interesting 
principal scale-independent differences (Amorim et al. 1997, Riisager et al. 2000). We 
will carefully compare relevant dimensionless parameters for both fields. An important 
dimensionless quantity is the density n multiplied by the third power of the scattering 
length a, i.e. na 3 , which is related to the classical mean free path A by na 3 ~ a/4irA. 
In general the size of A indicates which structure is preferred, e.g. mean field structure 
when A >• L and strong correlations when A ~ L (Bohr & Mottelson 1975, Siemens 
& Jensen 1987). This corresponds roughly respectively to na 3 <C 1 and na 3 ~ 0.1 — 10 
for atoms where cl/Lq ~ 10 — 100 and nuclei where a/Lo ~ 1 (away from resonance). 

The nucleon-nucleon interaction range in finite nuclei is comparable to the nuclear 
radius and the mean free path at low excitation energy (Bohr & Mottelson 1975, 
Siemens & Jensen 1987). These nucleon-nucleon properties arc not (necessarily) the 
same as in free space. The nuclear system is dense, and more than s-waves contribute. 
Both mean-field (single-particle or quasi-particle motion) and collective (rotational 
and vibrational motion of spatial, spin, and pairing) degrees of freedom are important 
in realistic descriptions (Bohr & Mottelson 1975). For atoms and molecules it is 
intuitively clear that the important degrees of freedom must be those related to the 
spatial orientation, the intrinsic excitations and the relative motion of the individual 
atoms (Haken & Wolf 1995). For both nuclei and atoms or molecules macroscopic as 
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well as microscopic properties are therefore necessary ingredients in the descriptions. 

The self-bound, leptodermous (thin skin) nuclear systems (Myers & Swiatecki 
1969) emphasize the importance of surface properties as for mesoscopic condensed 
matter systems (Reimann & Manninen 2002). Similar deviations from homogeneous 
bulk properties appear as correlations in trap- or optical lattice confined atomic 
systems (Bloch et al. 2008, Giorgini et al. 2008). The common theoretical quantitative 
techniques include the mean-field and BCS pairing approximations, linear response 
(including random phase approximations), ab initio variational minimization, shell 
model diagonalization in restricted Hilbert spaces, and density functional methods. 

The present paper is devoted to a comparison of common concepts and structures 
in nuclei and cold atomic gases. Striking similarities are seen in many discussions that 
employ the notion of condensation and paired states in both subfields. However, 
a careful comparison is demanded when concepts that traditionally refer to infinite 
systems are applied to finite systems such as nuclei or trapped atomic gases containing 
only very few atoms. Our aim is to facilitate and clarify exchange of important ideas, 
basic concepts and essential techniques between the two broad fields of physics. For 
a comprehensive review of experimental and theoretical developments in cold atomic 
gases we refer to (Bloch et al. 2008, Giorgini et al. 2008). 

Our starting point will be mean-field theory and we define the concepts of 
condensate and paired state in this context. Particular care is taken to address 
the problems with finite particle number and we therefore discuss at length the 
notion of condensate and condensate fraction. Paired states of fermionic particles 
in both nuclear and cold atomic gases are discussed mainly within the BCS theory 
which is widely used in both arenas. We also discuss the interesting physics of the 
transition from weakly-coupled paired states to a condensate of strongly bound two- 
body molecules as it is described within the BCS formulation (BCS-BEC crossover), 
and we address strongly- interacting Fermi gases in the unitarity limit. Along the way 
we take care to distinguish between generic features and model-dependent details that 
originate from certain convenient choices of model-space, effective interactions, etc. 

An important topic is the comparison of different scales in cold atomic gases and in 
nuclei. In particular, the isolation of relevant dimensionless quantities is essential. We 
undertake such a discussion next, explaining the degrees of freedom of both fields, and 
pointing out similarities and differences. An extended discussion of the cold atomic 
two-component fermionic gas is given from a nuclear physics perspective within the 
traditional shell model picture that should be helpful for applying nuclear physics 
techniques and concepts directly to cold atomic gases. As an example, we discuss 
few-body fermionic trapped atomic gases and compare them to light nuclei. 

We then turn to some specific examples of bosonic structures and condensates in 
few-body nuclear systems and also in trapped few-body atomic gases that have been 
studied theoretically in recent years. We employ a concept of quantum localization 
which can predict zero-temperature phases of quantum matter, and show how it is 
connected to more common measures of quantum degeneracy The case of finite 
systems of bosons is subsequently discussed with emphasis on the competition between 
localization of particles and mean-field states. For finite systems localization is 
unfavourable to the formation of a condensate. In nuclear physics, we apply these 
concepts to light nuclei that are believed to have cluster structures, and for which 
condensation of a-particles have been proposed as a potential structure of the system. 

Efimov physics is discussed next. Originating in nuclear physics, the concepts 
of an infinite number of bound three-body states around the threshold for two-body 
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binding which has recently been experimentally studied with cold atomic gases. This 
constitutes a beautiful example of crossdisciplinary fertilization. We discuss the 
concept of universal three-body bound states and give criterions for which systems 
are best suited for its observation. Finally, we discuss neutron-rich matter within the 
context of the BCS-BEC crossover. Again, we focus on few-body aspects and consider 
light halo nuclei that have a low-density surface occupied mostly by neutrons. Here 
we address which features related to the crossover can be described as universal and 
which are model-dependent. 

The current presentation is organized in sections with slightly different aims. The 
discussion of condensates in Section[2]is directed from atomic towards nuclear physics, 
implying that the content is well known in atomic and has to be defined properly 
in nuclear physics. The properties of pair correlations in Section [3] are common to 
nuclear and atomic physics but there are subtle issues with choices of model-space and 
residual pairing interactions in nuclear physics that we point out. Section [4] contains 
a description of atomic gas properties such as length and energy scales, degrees of 
freedom, etc., translated into nuclear physics language. Section [5] discusses few-body 
applications of the concepts of atomic boson gases in nuclear physics and in Section 6 
we specialize to the topical Efimov effect in both nuclei and atomic gases. Section 7 
is devoted to a discussion of BCS-BEC crossover physics in neutron-rich matter with 
particular emphasis on light nuclei and particular choices in the models that describe 
the nuclear interactions. Finally, section [8] contain a brief summary and conclusions. 



2. Condensates 



Bose-Einstein condensates were initially predicted as the result of a phase transition 
of non-interacting bosons below a critical temperature. This was at first considered 
a pathology of no interaction, but after the discovery of superfluidity in bosonic 4 He 
in 1938 by Kapitza (Kapitza 1938), Allen and Misener (Allen & Misener 1938), Fritz 
London had the insight that superfluidity might in fact be a sign of condensation 
(London 1938). A few years later Landau explained superfluidity in terms of the 
linear spectrum of the low-energy excitations in the system (Landau 1941). He also 
introduced a two-fluid model with a superfluid and a normal component, the former 
of which is related to the condensed particles. Superfluidity is connected with the 
excitation spectrum, whereas condensation is a property of the 'natural' stationary 
state of the system as given by the one- or two-body density matrix introduced below. 
From the point of view of atomic and nuclear physics, the latter concept is the more 
important and hence our focus in the following. 

As a large number of identical bosons are allowed simultaneously in the same 
single-particle quantum state it can be macroscopically occupied. An important aspect 
of condensed matter physics and the theory of phase transitions in three-dimensional 
space is that of long-range order. As we will discuss below, this concept is also useful in 
the case of finite systems where it is connected to the coherent state that a Bose system 
can condense into. The properties of boson condensates are extensively investigated 
in laboratories under a variety of conditions defined by ingenious combinations of 
external electromagnetic fields (Pethick & Smith 2008, Bloch et al. 2008, Grynberg 
& Robilliard 1999). In contrast, for fermions only one particle is allowed in a given 
quantum state. However, an entity consisting of an even number of fermions exhibits 
boson symmetry properties, and in particular these composite particles can occupy 
the same quantum state. The Cooper pairs of fermions in the BCS theory are 
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examples arising due to an attractive two-body interaction (Cooper 1956, Bardeen 
et al. 1957, Leggett 2006). Superconductivity and superfluidity are in this sense 
the same phenomena which can be seen in Landau's theory of phase transitions 
for these systems (Landau & Lifshitz 1958, Pitaevskii & Lifshitz 1980). For a 
summary of the concepts and historical development in superfluidity we refer to 
Leggett (Leggett 1999). 

In nuclei, prominent examples of bosonic substructures are a-particles which from 
the early days were suggested as building blocks of nuclei (Wigner 1937, Wheeler 1937, 
Wefelmeier 1937, Brink 1966), and pairs of nucleons within the BCS theory (Bohr 
et al. 1958, Bohr & Mottelson 1975). Numerous examples exist of both fundamentally 
bosonic atoms and molecules and also of binary bound fermionic molecules (Hakcn 
& Wolf 1995). When composite boson structures begin to get close to each other 
and interact, the ideal boson structure is no longer maintained. The intrinsic degrees 
of freedom become important and true microscopic fermionic nature is revealed. We 
shall here first discuss the definition and structure of boson condensates and in the 
next subsection the connection to bosons formed by two fermions via BCS-pairing. 

2.1. Structure of Bose condensates 

The ideal Bose condensate (BEC) is defined through quantum statistics for bosonic 
particles (Fetter k Walecka 1971, Pcthick & Smith 2008, Pitaevskii & Stringari 2003). 
The prediction is that below a certain temperature these systems will have a 
macroscopic number of particles, N , in the lowest energy state. Here macroscopic 
occupation refers to the thermodynamic limit where limjv->.oo {Nq/N) has to be finite 
when N is the total number of particles. The ideal Bose condensate with N = N 
is then by definition non-interacting since the inevitable correlations otherwise would 
prevent some particles from being in the same quantum state. This depletion of the 
condensate due to interactions depends on na 3 where a > is the scattering length 
characterizing the two-body interaction at low energy and n is the particle density. 
When a < 0, the bosons attract and only a limited number of particle can condense 
without collapsing as we will discusse below. However, to achieve condensation in 
an experiment the particles must be confined inside a trap. The presence of such 
an outer one-body trapping potential alters the particle spectrum but BEC can still 
occur for dilute gases in spatially extended traps (Pethick & Smith 2008, Pitaevskii 
& Stringari 2003). 

To be accurate in discussions of BEC we need to decide on an appropriate and 
robust definition. The most widely accepted basic definition, revolving around the 
density matrix for the system in question, was given early on in (Penrose 1951, Landau 
& Lifshitz 1958, Penrose & Onsager 1956, Yang 1962). We shall follow (Yang 1962) 
and specify properties of the TV-body density matrix arising from the wave function \I/ 
for the many-body system of N particles. We first define the one-body density matrix, 
Pi, by 



where r,,i = 1,..,N, are the particle coordinates. The choice of integration omitting 
i = 1 is arbitrary as \I/ is symmetric under any particle interchange. We assume that 
^ is normalized to N after integration over all coordinates of \^\ 2 . 




Pl (r,r<) = Tr\y(r))(y(r<)\ 
/ dr 2 dr 3 ...dr N y*(r,r2,...,r N )y(f*,r 2 ,...,r N ) , 



(1) 
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The positive definite matrix p\ with two continuous indices is hermitian and can 
thus be diagonalized by a unitary transformation, i.e. we can decompose this density 
matrix as 

i 

where are the eigenvalues (which are positive) and fa the corresponding single- 
particle eigenfunctions each normalized to unity. We order the eigenvalues to decrease 
in size with i. Integrating the diagonal part of the density matrix we get N, i.e. 

[ dr Px {r,?)=Y,h = N . (3) 

J i 

We define the system as a simple BEC when one and only one of the eigenvalues 
Ai remains of the order of magnitude A~ for increasing N, i.e. one large eigenvalue 
in the thermodynamic limit. Thus the wave functions for a BEC is dominated by a 
wave function that is a simple (symmetric) product of the same single-particle wave 
function fa, i.e. ^ s» Tiifa(ri). Notice that this holds for condensates that are close 
to the ideal non-interacting case. A counterexample is atomic 4 He which has very 
strong correlations and thus only 10% of a given sample will be in the condensed 
state. This definition inherently compares structures for different (large) numbers of 
particles. More than one eigenvalue may be large even in the limit of large N. Such 
states are called fragmented BEC and they have recently received increasing attention, 
particularly in connection with gases of more than one species of atoms and also with 
respect to condensation of asymmetric (different number of particles) two-component 
systems (Pethick & Smith 2008, Pitaevskii & Stringari 2003, Mueller et al. 2006). 

The concept of off-diagonal long-range order is closely related to the behavior of 
the one-body density matrix in the limit where \r — r"| becomes large compared to the 
system size (Yang 1962). Let us then consider a possible limiting form of pi, i.e. 

lim pi = F*(r)F(f*)+p' 1 (r,r') , (4) 

| r — f' | — > oo 

where the first term has the product form obtained by diagonalization. If the 
function p' x tends to zero in this limit one would say that the system is a BEC with 
order parameter F which depends on the coordinate. Then F is interpreted as the 
condensate wave function, F — > faVN, containing A" particles and obtained in the 
mean-field approximation. Thus Eq. Q separates p\ into the product part related to 
the largest eigenvalue and the rest of the sum in Eq. ^ denoted by p[ . One can show 
under rather general conditions that p' x vanishes for homogeneous 3-dimensional Bose 
systems (Leggett 2006) and also in the case of trapped condensates (Naraschewski & 
Glauber 2001). Thus this alternative definition of BEC is equivalent to that of one 
large eigenvalue in the expansion in Eq. 

The one-body density matrix for a BEC in a trap was experimentally measured 
about a decade ago (Bloch et al. 2000) and the expected coherence was confirmed. This 
corresponds to what is known as first-order coherence (Naraschewski & Glauber 2001). 
Second-order coherence (related to the two-body density matrix) is connected to 
the famous Hanbury-Brown and Twiss effect (Hanbury-Brown & Twiss 1956, Baym 
1969, Baym 1998). The corresponding bunching effect has been observed in optics 
(Bachor & Ralph 2004) and in cold bosonic atoms (Yasuda & Shimizu 1996, Foiling 
et al. 2005, Schellekens et al. 2005). The analogous anti-bunching in electronic systems 
(Henny et al. 1999, Oliver et al. 1999, Kiesel et al. 2002), for neutrons (Iannuzzi 



CONTENTS 



8 



et al. 2006), and in fermionic cold atoms (Rom et al. 2006). Recent experiments 
have even begun to study third-order coherence (large eigenvalues in the third-order 
correlation function) and the nucleation of a 4 He BEC (Hodgman et al. 2012). 

2.2. Center- of -mass motion 

There is one possible flaw in this equivalence between the two definitions of BEC 
above. Unfortunately that is related to trapping in external one-body potentials 
which is employed in all laboratory experiments. The problem is that the center- 
of-mass motion is not decoupled as for a self-supported system. The center-of-mass of 
the N particles do not coincide with the center of origin for the external field. Then the 
particles can for example move coherently in rotational states without changing their 
relative motion. This would show up as simultaneously existing BEC-like structures 
with eigenvalues of comparable magnitude even when the number of particles increases. 
This was pointed out in (Pethick & Pitaevskii 2000) and an alternative definition of 
BEC was suggested to remove the effect of the center-of-mass motion, i.e. using 
the internal one-body density matrix where only relative coordinates enter. Then 
different rotational motion build on the same intrinsic BEC would all lead to one 
large eigenvalue of the internal one-body density matrix. 



10 12 14 16 18 

N 



Figure 1. Condensate fraction, ct, as function of the number of particles, N, for 
a gaussian mean-field wave function when the center-of-mass has been separated 
out as discussed in the text. 

We can illustrate by a gaussian mean-field wave function of range b where 
the center-of-mass coordinate separates out leaving only relative coordinates in the 
internal wave function, i.e. 

* in *({r,}) = (&^r 3(Ar ~ 1)/2 exp( V/(26 2 )) , (5) 

N 1 W 

p 2 =J2<£ > n = n-R , R= ^5>* ' ( 6 ) 

»=1 i=l 

where the coordinates qi are measured from the common center-of-mass R. This wave 
function is invariant under rotations around R. Following (Pethick & Pitaevskii 2000) 
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the internal one-body density matrix p(q 1 ,q[) is now found by inserting qjv 
-Eili^i in Eq. (i}, i.e. 



p(Qi,Q'i)^ J d 3 q 2 d 3 q 3 ...d 3 q N _ 1 \ 



( Q 2 + q' 2 , {N-2){q' + qy 

The condensate fraction obtained through the largest eigenvalue is then (Gajda 2006) 
Cf = 8/(1 + - 2/N) 3, which decreases with N from 1 for N = 2 towards about 0.57 
for large N. This is illustrated in figure [I] However, the choice of relative coordinates 
is arbitrary (Gajda 2006) and we could as well choose gi supplemented by a set 
of N — 1 independent Jacobi coordinates. Then the density matrix corresponding to 
Eq. ^ would factorize and give Cf = 1. This fact has also been addressed in (Yamada 
et al. 2008, Yamada et al. 2009). Clearly these discussions are only relevant for finite 
systems where the center-of-mass can move. Thus the distinctions become less and 
less important with increasing particle number, but at the same time interesting for 
a smaller number of particles in BEC-like states. 



2.3. Mean- field description 

The initial definition of BEC as a number of particles in the same single-particle 
state immediately points to a mean-field approximation or independent particle model. 
The only assumption is that the many-body wave function is a symmetric product of 
single-particle wave functions. With this restriction the variational principle gives the 
lowest energy solution for any Hamiltonian. This implies that the average effect of 
the interactions between the particles are included in this approximation. Obviously 
any interaction would prefer to correlate the particles either to avoid or to exploit the 
interaction. This cannot be expressed in the mean-field product wave function which 
is the ideal BEC structure. Thus any deviations beyond the mean-field must reduce 
the BEC content. In other words interactions introduce correlations which attempt 
to destroy the BEC structure. 

For cold atoms the zero-range approximation for bosons has been used 
successfully in the form of the self-consistent mean-field Gross-Pit aevskii equation 
(Gross 1961, Pitaevskii 1961, Pethick & Smith 2008, Pitaevskii & Stringari 2003). 
It is demonstrated for dilute systems to be very efficient in reproducing data for 
interacting condensates with inter particle repulsion. For attractive interactions 
(a < 0), the system collapses unless some measures are taken to avoid this part 
of the Hilbert space. However, quasi-stable solutions exist of finite size and energy if 
N\a\/b < 0.65 where b is the trap length of the external harmonic potential (Ruprecht 
et al. 1995, Kagan et al. 1998, Pitaevskii & Stringari 2003, Pethick & Smith 2008, Bohn 
et al. 1998, S0rensen et al. 2004). Figure [2] shows the effective potential barrier of 
an attractively interacting condensate based on a gaussian variational calculation 
including also higher-order effective range corrections (Zinner & Th0gersen 2009). 
The variational wave function is 

where b — \Jh/mu> is the oscillator length, N is the particle number, and q is the 
dimensionless variational parameter. The quasi-stable solution are located in the local 
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minimum around q <~ 1. These solutions decay into deeper lying (collapsed) states 
of the many-body system, e.g. non-condensate like states where a number of binary 
or many-body bound states are present perhaps in some mixture with dilute atomic 
states confined by the external field. When N\a\/b is much lower than the critical 
value the lifetimes of the condensate is much longer than the experimental time. 




Figure 2. Energy of a BEC with fixed N\a\/b as function of the variational 
parameter q, i.e. the size of the condensate (proportional to b). The higher-order 
interaction term g2 = a 3 /3 — ar e /2 includes the effective-range correction through 
r e . It modifies the height and shape of the barrier through which macroscopic 
quantum tunneling occurs. The left panel has N\a\/b = 0.5 < 0.65 and the right 
panel shows results for N\a\/b = 0.7 > 0.65, i.e. on both sides of the stable 
regime when g2 can be neglected. The right panel demonstrates that higher- 
order corrections can stabilize an attractively interacting condensate. Taken from 
(Zinner & Th0gersen 2009). 

The mean-field treatment provides the optimal solution where the BEC product 
wave function is maintained. The most sensible definition of correlations in many-body 
physics is in terms of deviations from the simple mean-field product wave function. 
Therefore when we define BEC through the one-body density and its eigenvalues 
we have an intuitive access to understand the effects of correlations on a BEC. 
Expanding the full correlated many-body wave function on a basis of different mean- 
field product wave functions, a number of non-zero eigenvalues (A^) would appear with 
corresponding single-particle wave functions. The condensate fraction is then given by 
the largest eigenvalue, A*, as A* /N . Imagine that we turn on interactions slowly in an 
ideal BEC. The condensate fraction then decreases from unity and many eigenvalues 
of order 1 would appear. Approaching the thermodynamic limit by increasing N while 
maintaining the same interactions would again increase the largest eigenvalue and the 
condensate fraction towards unity. Even rather strong two-body interactions would 
not destroy the BEC structure reached for large N (Mueller et al. 2006, Th0gersen 
et al. 2007). This does not necessarily mean that the BEC structure is unaltered by 
the interactions and the additional particles. It means that the best approximation 
by a product wave function, depending on interactions and N, is approached (Licb 
et al. 2002). 
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2-4- Condensates of Pairs of Fermions 

Combining two fermions into a, possibly bound, entity produces a system with boson 
characteristics. Redefining the degrees of freedom from single-particle to two-particle 
center-of-mass and relative coordinates, and integration over the relative coordinates in 
the density matrix, allow use of the general definition of boson condensates. However, 
if the intrinsic two-body structure cannot be frozen, we need a generalization of the 
BEC concept to systems of fermions. Due to the Pauli exclusion principle any one- 
body density matrix will always have eigenvalues that are less than or equal to one, 
so using this to look for coherence makes little sense. However, as Yang (Yang 1962) 
shows, the two-body density matrix, p 2 , for fermions should be used instead of the 
one-body density matrix p\. The definition of p 2 is 

P2(n,f 2 ,r' 1 ,r' 2 ) (9) 



dr 3 ...dr N ^ (f u r 2l ...,f N )^f(f' 1 ,f 2 , ...,r N ) , 

which may have eigenvalues, A^, of order N, i.e. 

P2(ri,f 2 ,r' 1 ,f^) =2^Ai0*(ri,f 2 )0i(r'i,f|) . (10) 

It is then natural to find these eigenvalues and define a fermion condensate as the 
structure where one of these eigenvalues is large as for BEC. The connection to a 
boson condensate of composite particles is then readily seen by change to relative and 
center-of-mass coordinates r = f\ — f 2 , R = (f\ +r 2 )/2, and subsequent integration 
over the "intrinsic" r-coordinates, i.e. 

pi{R,R') = J p 2 {r 1 ,r 2 ,r 1 ,r 2 )drdr . (11) 

If (f>i(r\, ity factorizes in the (R,r) coordinates the fermion condensate reduces to a 
boson condensate of pairs of fermions. This two-body density matrix determines the 
second-order coherence of the system for both fermions and bosons as was mentioned 
above. 

In many cold atomic gas experiments with fermions the system has two internal 
states, similar to electrons in condensed-matter systems. The fermions carry an 
effective half integer quantum number (the physical origin of this will be discussed in 
later sections). We have suppressed this index here. When internal states are present, 
one can think of the coordinates in Eq. [9] as containing both spatial coordinate and 
internal state information. 

Some very simple, although pathological, limits of very few particles have 
been discussed in nuclei. For two identical bosons the one-body density matrix is 
automatically of the form in Eq. ^ with only one term in the summation which means 
maximum eigenvalue and ideal Bose condensate. Thus when the center-of-mass wave 
function is controlled by an external field or included as in (Pethick & Pitaevskii 2000) 
this definition implies that two bosons always form a condensate. Two fermions 
with identical orbital wave function and antisymmetric relative spin wave function 



automatically factorize as in Eq. ( 10 1 with one term. This implies that two identical 
fermions in a state of lowest energy are in a fermion condensate. Furthermore, the 
composite system of two fermions is by definition in a boson condensate consisting 
of one particle. These rather pathological limits demonstrate that composites of few 
particles with properly chosen definitions can fulfil the criteria for a condensate. Of 
course one must still determine whether this has any additional physical significance 
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in the relevant system at hand. We will discuss these questions in a few simple cases 
in the second part of the review. 

3. Pair correlations 

Before we discuss pairing correlations in two-component Fermi systems, we first 
address an important question related to the choice of model space and interactions. 
This is necessary since there are several choices involved before one arrives at some 
properly defined model to which pair correlations and a particular method of solving 
pairing Hamiltonians can be applied. This is the case for both finite and infinite 
systems, and it is important to keep track of the effective interactions and model 
spaces that are applied when addressing generic features or universal behavior. 

In the spirit of density functional computations, a sufficiently well chosen 
density functional would provide any desired accuracy of the true energy (Kohn & 
Sham 1965, Dreizler & Gross 1995). However, correlations are entirely missing in 
the one-body structures of the wave functions. To investigate effects of correlations 
it is therefore necessary to use a Hilbert space that goes beyond one-body product 
states. As correlations are driven by interactions it is also necessary to know precisely 
which interaction to use in the extended space. These connections and the practical 
implementations in nuclear and atomic physics are completely different. 

A natural starting point in both subfields is some form of self-consistent mean- 
field model. For fermions the most common choice is the Hartree-Fock approximation 
(Fetter & Walecka 1971, Negele 1982, Siemens & Jensen 1987, Pethick & Smith 2008). 
For identical bosons (fermions) the Hilbert space is (anti)symmetrized products of 
single-particle wave functions. Several extensions of the Hartree-Fock method are 
in use. In nuclear physics, a prominent example is the Hartree-Fock-Brueckner 
scheme (Fetter & Walecka 1971, Siemens & Jensen 1987) which takes medium effects 
of Pauli blocking into account in a self-consistent manner. For cold atomic gases 
similar medium effects have been taken into account in various other ways (Heiselberg 
et al. 2000, Stoof et al. 2009). For the present discussion we keep things simple and 
stay at first within the framework of the standard Hartree-Fock approximation. Later 
on we comment on some important features of going beyond this approximation in 
connection with polarization effects, induced interactions and the renormalized zero- 
range BCS model. 

A mean-field Hartree-Fock calculation requires suitable two-body interactions 
as input. In atomic gases the diluteness of the system makes the zero-range 
approximation extremely accurate even around broad Feshbach resonances where 
the effective range remains small (Kohler et al. 2006) (a more detailed discussion 



of Feshbach resonances is given in Section 4.1 ). However, it is remarkable that nuclear 



physics mean-field calculations employ zero-range interactions since nuclei are not 
dilute as the range of the nucleon-nucleon potential is comparable to the interparticlc 
distance in the nucleus. Cold atomic gases are usually dilute in this sense and zero- 
range interactions are more appropriate (at least away from interaction resonances as 
we will discuss below). An interesting comparison can be made to fermionic liquid 
3 He which is a dense system. 

In general, it is crucial to relate the chosen interaction with the allocated Hilbert 
space (Barrett 1980, Suzuki & Lee 1980, Suzuki 1982, Brandow 1967, Poves & 
Zuker 1981). This must be done whether one uses zero-range or some other form 
as the model interaction, although using a finite-range potential can often help avoid 
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divergences. The physical interaction must be transformed to apply in a restricted 
Hilbert space, and the operators for other observables than the energy must be 
correspondingly transformed. Properly done all physics is then correctly maintained. 
The connection between basic and effective (transformed) interactions can easily be 
complicated or downright impossible to trace. It is for example rigorously necessary to 
work with up to TV-body interactions, not only the simple initial two-body interactions 
(Suzuki & Lee 1980, Suzuki 1982). In any case, nuclear and atomic physics differ 
tremendously in both effective interactions and Hilbert space even when identical 
techniques are employed (Barrett 1980). 

The nuclear force is at low energy expanded in terms of the relative momentum 
between the nucleons, and only terms up to second order is maintained (Siemens 
& Jensen 1987). This Skyrme interaction corresponds to zero-range in coordinate 
space (Skyrme 1956) . The inevitable collapse for attractive potentials in Hartree-Fock 
applications is avoided by a completely phenomenological (still zero-range) density 
dependence, or alternatively multi-body forces (Negele & Vautherin 1972, Negele & 
Vautherin 1975). The strengths of the different terms are then adjusted to reproduce 
observables like energies and sizes of a series of nuclei. The connection to the basic 
nucleon-nucleon interaction is lost. 

In atomic physics the zero-range interaction in mean-field calculations leads for 
identical bosons to the Gross-Pitaevskii equation (Gross 1961, Pitaevskii 1961, Pethick 
& Smith 2008, Pitaevskii & Stringari 2003). The total wavefunction for N bosons in 
factorized form, ^> = Tlf =l <t>(ri) , leads to the stationary Gross-Pitaevskii equation for 
<f> which reads 

- ^^(r) + Wr>(r) + ^l^l V(0 = E^f), (12) 

where m is the boson mass and E the energy per particle. The non-linear interaction 
term, which gives rise to many interesting physical effects in condensates, are well- 
described by this equation (Pitaevskii & Stringari 2003). The strength of the 
interaction is chosen to reproduce large distance, or equivalently low-energy scattering 
properties within a confining external field. This amounts to reproduce the physical 
atom-atom scattering length by the Born approximation of the potential (Pethick & 
Smith 2008, Pitaevskii & Stringari 2003). The intuitive explanation is that Hartree- 
Fock product wave functions in the one-body external potential corresponds to free 
uncorrelated solutions for which the Born approximation is sufficient. For repulsion 
the solution can then immediately be obtained. The collapse for attraction is avoided 
by restriction of the Hilbert space to larger distances or, equivalently to, lower energies. 
The connection to low-energy two-body scattering properties is maintained, and the 
strength related to the scattering length in this limit. 

3.1. The BCS-approximation 

Correlations in general and pair correlations in particular arise only by going 
beyond the Hartree-Fock or Gross-Pitaevskii approximation. The neglected residual 
interaction must be defined relative to the main part accounted for in the mean- 
field approximation. In atoms the residual interaction is rather well-determined as 
connected to the basic physical interaction. However, for nuclei this immediately 
presents the problem that the procedure is not unique since the starting point is 
phenomenological. The residual and mean-field interactions for nuclei must then be 
related through the same phenomenology. Conclusions about physics beyond the 



CONTENTS 



14 



region where the parameters are adjusted can be very uncertain. Recent advances in 
deriving low-momentum interactions for use in few- and many-body nuclei promise 
to improve on this situation by decoupling the troublesome high-momentum parts of 
the typical nucleon-nucleon interaction while preserving the correct scattering phase 
shifts (Epelbaum et al. 2009, Bogner et al. 2010). These developments may be relevant 
for the interchange between nuclei and cold atoms. However, as these technical 
issues are not necessary for the purpose of discussing pairing correlations in the BCS 
approximation, we focus here on the more phenomenological approach traditionally 
employed for nuclei. 

Pair correlations can be incorporated on an equal foot with the single-particle 
degrees of freedom by extension of the Hartree-Fock to the Hartree-Fock-Bogoliubov 
(or Bogoliubov-de Gennes in the atomic literature) approximation (de Gennes 1999, 
Bohr & Mottelson 1975, Ring & Schuck 1980, Dobaczewski & Nazarewicz 2012). To 
grasp the main idea, the simplest extension is sufficient which is to supplement the 
mean-field treatment by the BCS approximation. For simplicity we assume identical 
energy spectra, e, for the different internal states of the fermions we have in mind. In 
traditional superconductivity, the two internal states of the electrons are of course the 
spin degrees of freedom, whereas in a two-component ultracold atomic Fermi gas the 
internal states are typically two states with different hypcrfinc projections (this will 
be discussed in greater detail below). We will also assume that the same number of 
particles occupy each internal state forming pairs. The latter condition means that 
we consider here only such balanced systems where the original implementation of the 
BCS theory applies. When the system is imbalanced, a number of exotic states like 
the FFLO pairing state (characterized by a oscillating pairing gap parameter) have 
been predicted (see (Casalbuoni & Nardulli 2004) for a recent review). In nuclei, the 
different internal states are given by time-reversal symmetry (Bohr et al. 1958, Bohr 
& Mottelson 1975, Ring & Schuck 1980, Siemens & Jensen 1987), and in cold atomic 
gases by the two decoupled hyperfine states with controlled frozen occupation (see 
Section [4| . The role of time-reversal symmetry was already discussed in the context 
of superconductors by Anderson (Anderson 1959). The Hamiltonian, after subtraction 
of the Lagrange multiplier term /iN, is 



where the single-particle energies e i7 the creation <zj, and annihilation a,j operators 
refer to the single-particle mean-field states \i). The internal states, time-reversed or 
hyperfine, are denoted by a bar and we assume Ci = ej. The single-particle energies are 
measured with respect to the chemical potential /i multiplying the number operator N. 
If the spectra of each component are not the same, we need a sum of two terms each 
related to the two types of internal states. The matrix elements Gij of the residual 
interaction Hr, beyond the mean-field H m f then denotes the two-body interaction 
between pairs of particles in time-reversed (hyperfine) mean-field states. This piece 
of the interaction is called the pairing interaction, and it should not be confused with 
neither the full two-body interaction employed in the mean-field potential nor any 
other residual contribution to the matrix elements, e.g. related to other correlations 
like vibrational or rotational excitations in the RPA treatment. 



Hp = H mf — fj,N + Hr 



(13) 
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The BCS approximation then consists in 

fr„«fl m /-/*iV-5^A 4 (oJoi + a 5 Oi) , (14) 

i 

which is diagonalized exactly by the Bogoliubov transformation 

a] = ma] - vioq , at = Uiat + UjO; , (15) 

where the new quasi-particle creation a\ and annihilation on operators also obey 
fermion anti commutation rules. The result is 

w U + E i( a Ui + a \ a *) > ( 16 ) 

i 

i i 

The occupation numbers vf and quasi-particle excitation energies Ei are 

«?=!-«?= 2^), (18) 



E, 



^/(e t - M )2 + A? . (19) 



The chemical potential /i and the gaps Aj are determined for a given average number 
of particles Nq = (N) from the equations 

A* 



AT = 2W , 2A^Vg^. (20) 



The new BCS ground state wave function related to the aj quasi-particle operators is 
expressed by 

\BCS >= Yl % { Ul + «<ajoj)]0 > , (21) 

where |0 > is the "vacuum" related to the a% operators. It is also possible to write an 
analogue of the BCS theory presented here that conserves particle number by using 
projection (Dietrich et al. 1964, Adornes & Kyotoku 1992, Leggett 2006) although 
these can become considerably more complicated. Schematic and exactly solvable 
pairing models have also been studied in great detail (Dukelsky et al. 2004, Pan & 
Draayer 1999). 

The pair correlations obtained in the BCS-approximation are derived from the 
residual interaction expressed in terms of the two-body matrix elements, Gy, of the 



interaction between pairs in time-reversed states, see Eq. (13) An approximation 
of constant pairing matrix-elements, G, independent of i and j is very instructive 
and often also reflects correct average pairing properties. This is an approximation 
when the single-particle states are different from plane waves even for a zero-range 
interaction. Then the gap Aj = A also becomes state-independent and the occupation 



v~ in Eq. ( 18 ) are seen only to deviate from zero or unity around the Fermi 
surface where the least bound particles reside. This indicates a tendency of pairing as 
a surface phenomenon for finite systems (Bohr et al. 1958, Bohr & Mottelson 1975). 

A zero-range residual interaction leads to constant matrix-elements for plane wave 
single-particle states. However, the summations for constant G, e.g. in Eq. |20] ), 
produce an unphysical divergence which is removed by a cut-off in the sums at a finite 
value of i. The price is an adjustment, depending on the cut-off, of the value of G to 
a physical observable (Huang 1987, Pethick & Smith 2008). The mean-field solutions 
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usually deviate from plane waves arising in free space with no external confinement. 
However, the gap equation can in fact be solved for arbitrary mean-field states with 
constant matrix-elements, as discovered early on by Richardson (Richardson 1963). 
Any residual pairing interaction, including zero-range, will have state-dependent 
matrix-elements in general. It is then as convenient to use a realistic residual 
interaction of finite range which automatically removes the divergence. Numerical 
treatment is then necessary but nowadays only marginally more difficult. 

The new ground state in Eq. (21 ) describes a system of (Cooper) pairs of particles 
with opposite momenta and zero total angular momentum. This result applies to 
continuous mean-field spectra for any attractive short-range residual interaction, and 
for discrete (finite) systems above a critical strength, G c (Bohr et al. 1958, Fetter & 
Walecka 1971, Siemens & Jensen 1987). This ground state is sometimes called the 
pair condensate to emphasize the analogy with boson condensates. The first excited 
state is, for constant matrix elements, at 2A corresponding to the energy required to 
break one pair. The gain in energy from the normal state with a Slater determinant 
wave function (a filled Fermi sea) to the superfluid state is g(fi)A 2 /2 where g(fi) is 
the single-particle level density at the Fermi energy fi (Fetter & Walecka 1971). Here 
we ignore low-energy collective modes which appear at zero excitation energy in the 
infinite system limit but which are pushed up to finite energy (i.e. they have an 
excitation gap) in finite systems as well (Landau & Lifshitz 1958). 

These conclusions are quite remarkable when considering that in three dimensions 
it takes a finite attraction to even produce a bound state. Here (almost) any small 
attraction produces a qualitatively different solution. In the original treatment of 
the effect by Cooper (Cooper 1956) this occurs since the pairing is assumed to occur 
only near the Fermi surface. The problem then becomes effectively two-dimensional 
or, equivalently, a one-dimensional problem with a linear spectrum around the Fermi 
energy. Here it is well-known that an arbitrarily weak attraction leads to a bound state 
(Landau & Lifshitz 1981, Simon 1976, Volosniev, Fedorov, Jensen & Zinner 2011). The 
presence of the Fermi sea, however, implies that we are dealing with a true many-body 
effect. In the cold atomic gases the zero-range approximation is accurate and should be 
used also in the Cooper pair problem. However, the pairing interaction is now constant 
everywhere in the Fermi sea and this introduces the need for renormalization of the 
divergent bound-state equation. The bound-state remains but its energy changes. 
In nuclear physics, the problem is usually treated by employing a cut-off either by 
hand or by using a finite-range interaction that gives vanishing matrix elements above 
some energy scale. Thus, we see that different methods are in fact needed to extract 
the same physics in different fields. We also caution that one should be careful to 
distinguish between the Cooper pair energy and the BCS gap parameter which have 
similar expressions as function of the pair interaction strength (or scattering length) 
but different exponents (de Gennes 1999, Leggett 2006). All this should serve as a 
warning of the trouble that one faces when extending these concepts to small systems 
with few particles both in nuclear and cold atomic systems. 

To summarize for nuclei, both mean-field and residual interactions are often 
chosen phenomenologically. It is then of less importance whether parameterizations 
are for a pairing residual potential, for the matrix elements themselves, or if less 
details are required, for the constant matrix elements combined with a cut-off. For 
atoms the mean-field potential is to a large extent determined by the external one- 
body potential although the physical two-body interaction also may contribute. The 
contributing matrix elements in mean-field and residual pairing interaction are not the 
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same although for atoms both are obtained directly from the low-energy scattering 
properties of the physical two-body interaction. It should be emphasized once again 
that the two distinct states that have the paired structure are time-reversed orbits in 
nuclei but they are states of different hyperfine projection in the case of cold atoms 



(we give more details in Section 4.3 1 



For atoms the procedure for the BCS-calculations could be to choose a potential 
of for example gaussian shape and adjust range and strength to reproduce the desired 
scattering length for particles in paired states. The range is very small compared 
to the dimension of the systems and the simplifications of a zero-range interaction 
is therefore tempting. However, this immediately introduces a divergence in the 



gap equations (Eq. (20) with Gij = G), and a renormalization procedure becomes 
necessary. The method usually adopted is based on the observation that the scattering 
length expression from the Lippmann-Schwinger equation formally resembles the gap- 



equation, Eq. (201, with the same type of divergence (similar to the procedures used 
for Fermi gases with repulsive hardcore interaction (Fetter & Walecka 1971, Pitaevskii 
& Stringari 2003)). Subtraction of these two divergent equations replaces the pairing 
matrix element by the physical scattering length for particles in paired states, i.e. 

2rf a = ? £k "^' (22) 

k k 



which combined with the number equation in Eq. (20) provide gap, A, and chemical 
potential, /i (Eagles 1969, Leggett 1980, Engelbrecht et al. 1997, Papenbrock & 
Bcrtsch 1999, Chen et al. 2005a). This derivation is based on the assumption that the 
mean-field wave functions are plane waves leading to constant pairing matrix elements 
for the zero-range residual interaction. Thus the approximation is exact for a non- 
interacting system in an infinite volume. To be valid the system must therefore be 
dilute and spatially extended. If we start with a finite-range potential then the limit 
of zero-range must be approached after the infinite volume limit. 

In the discussion above we have ignored the effect on pairing from polarization of 
the background Fermi gas. This effect was considered shortly after the introduction of 
BCS theory by Gorkov and Melik-Barkhudarov (Gorkov & Melik-Barkhudarov 1961) 
(often refered to as the GMB correction) who calculated a reduction by a factor of 
(4e) -1 / 3 ~ 0.4514 on the transition temperature from paired to normal phases in a 
dilute Fermi gas. The zero-temperature value of the gap discussed above is reduced 
by the same factor. This effect was later discussed in the specific context of ultracold 
atomic Fermi gases (Heiselberg et al. 2000) and similar effects have been pointed 
out for Bose-Fermi and Bose-Bose mixtures (Pethick & Smith 2008). Such medium 
effects are relevant for many condensed-matter systems such as for instance liquid 
3 He where the equilibrium ground-state is determined from a manifold of degenerate 
possiblities by medium effects (Vollhardt & Wolfle 1990). In the context of nuclear 
physics, this discussion also has a long history which is intimately connected with 
Hartree-Fock-Bruckner theory and the independent pair approximation (Fetter & 
Walecka 1971, Siemens & Jensen 1987, Dickhoff & Neck 2005). For nuclear matter, 
the inclusion of medium polarization reduces the gap similar to the GMB corrections 
(Chen et al. 1993, Wambach et al. 1993, Schulze et al. 1996). 

3.2. Pertinent properties of the BCS solution 

The BCS solution reveals two-particle correlations through the so-called pair wave 
function ^ pa ir, he. the amplitude for finding one particle in f\ and another in r*2. We 
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have 

4W(ri,f 2 ) cx (BCSl^in^if^BCS) , (23) 

where ijy is the spatial creation operator which is expressed by the mean-field wave 
functions fa in the external field as 

i 

In terms of the u and v coefficients of Eq.([T5| (Leggett 2006), one obtains 

*pair(6.,^2) OC J~] UjVj (fa (fx ) <fo (f 2 ) + 4' (^2)^(^1)) , (25) 

i 

which is a symmetrized product wave function of the spatially identical wave functions 
fa and fa. We have omitted the antisymmetric internal wave function, spin singlet for 
nuclei and antisymmetrized product of hyperfine states for atoms. 

The uu-factor, UjU, = Aj/2-E,, implies that the states within an energy interval 
of A around the Fermi level are the main contributors to ^ pa ir- When the mean- field 
wave functions are plane waves as for a homogeneous system the pair wave function in 



Eq. (231 becomes the Fourier transform of the uu-product with respect to the relative 
coordinate, and as such it depends only on the difference between ft and r*2. For 
finite systems like nuclei the center-of-mass motion of the two particles are equally 
important in ^ pa ir as suggested by the relation 2(rf + ff) = (f\ — f 2 ) 2 + (ri + r 2 ) 2 . 
The wave function separates exactly into an independent center-of-mass factor for a 
harmonic oscillator potential, and approximately for other potentials. 

The pair wave function, ^ pa in carries information about the two-body correlation 
for particles around the Fermi level. It can be used to compute properties related to 
the pairs as if it is an ordinary wave function. For instance, the relative mean square 
radius is obtained by 

(r 2 > = / d 3 nd 3 f 2 (fi - f 2 ) 2 |# pQir (ri,r2)| 2 , (26) 



where the center-of-mass coordinate in ^ pa ir can be removed completely for a 
homogeneous system and decoupled for a harmonic oscillator external mean-field. 



The expectation value in Eq. ( 26 ) is often called the mean square radius of a Cooper 
pair (Leggett 2006). 

For a homogeneous system of fermions with two internal states and a zero-range 
two-body pairing potential, ^ / pa ir(^) is for large values of r — \fi — f^l given by 
(Leggett 2006) 

^(r)oc Si ^e-^' ! (27) 
kpr 

where kp is the wavenumber at the Fermi energy and £' = h 2 kp/(mA) with A as 
the pairing gap at the Fermi surface. Thus £' is seen as the characteristic length 
of the pair correlation (historically one considers £ = £'/tt which is the so-called 
Pippard coherence length (Fetter & Walecka 1971)). The constant of proportionality 



in Eq. (27) depends on A and the density of state at kp. The generic form of the 
pair wave function is shown in figure [3] From Eq. ( 27 ) we clearly see that the value of 



kp determines the density of nodes of the pair wave function. These oscillations are 
often used as a signature of BCS behavior. However, in finite systems they can arise 
from other mechanism as we shall see in the nuclear case in Section [7j Except for the 
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exponential factor, the expression in Eq. (27) is exactly that of a free particle (or pair 
in relative coordinates) having an energy corresponding to kp (also shown in figure 
|3| . It is therefore apparent that the true many-body character of the BCS solution is 
found in the exponential decay of the pair correlation. However, in nuclei £' is large 
and therefore the exponential region is not entered before finite-size effects set in. We 
will return to this point when we consider BCS in finite nuclei in Section [7] 




5 10 15 20 25 30 35 40 45 50 



Figure 3. Plot of the pair wave function, Eq. \27\ , multiplied by kpr as 
function of kpr (solid (blue) line) for kp f = V200 (weak-coupling BCS limit). 
For comparison, we also show the non-interacting wave function with momentum 
kp (dotted (black) line) and the bound state wave wavefunction corresponding to 
the exponential damping in Eq. \27\ (dahsed (red) line). 



For finite systems large £' complicates the signature of BCS from the pair- 
wave function since there will always be an exponential suppression of the wave 
function for radii that are outside the potential range. Reduction of the BCS many- 
body wave function to a two-body pair wave function yields Eq. ( 25 1 which has the 
largest amplitude when the uw-factor is maximum, i.e. around the Fermi level where 
u 2 = v 2 = 1/2. This largest component (k — kp) is responsible for the oscillating 
term in the approximation of ^ pa ir in Eq. (27). The exponential damping on the 
other hand becomes pronounced at large distances and grows with oc A. bmce 
the condensate fraction (of Cooper pairs) is also propertional to A (Leggett 2006), we 
conclude that this damping is caused by the presence of the condensate at low energy 
(corresponding to large distance). 



3.3. The BCS-BEC crossover 

The physics of the BCS-BEC crossover originates from considerations of strongly 
coupled superconductors by Eagles (Eagles 1969) and atomic Fermi gases by Leggett 
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(Leggett 1980). The central observation is that by varying the two-particle interactions 
between the electrons/alkali atoms from weak to strong, one should observe a 
transition from the BCS paired state to a molecular BEC state where the relative 
two-body bound state is the relevant degree of freedom, hence the name BCS-BEC 
crossover. Intermediate stages would presumably have both pairs, molecules, and 
single fermions. It is important to emphasize that this is not a phase transition 
in the usual sense of Landau's theory of symmetry breaking transition but rather 
a smooth transition. The BCS-BEC crossover is the focus of intensive studies 
currently, both in condensed-matter, atomic, and nuclear physics. In particular, 
the use of optical lattices in ultracold atomic gas physics hold great promise for 
achieving insights into models that are typically applied in condensed-matter systems 
and are likely a key to a better understanding of high-T c superconductivity (Chen 
et al. 2005 a, Leggett 2006, Bloch et al. 2008, Giorgini et al. 2008). 

The physics associated with crossover, i.e. the transition from extended pairs 
to tightly-bound molecules, is found already in the relatively simple Leggett model 
(Leggett 2006) , which is expressed in the homogeneous version of the gap and number 
equations in Eqs. (20) and (22). These equations yield gap, A, and chemical potential, 
/i, for any given density n — k F /(3ir 2 ), where kp is the Fermi momentum related 
to the Fermi energy ep = h 2 k F /(2m). By dimensional analysis, A = epf(kpa) 
and /i = epg(kpa), where / and g can be determined numerically. We also have 
(kpa) 3, oc na 3 (typical values are given in Table (Jl|). In the ultracold gases one 
fixes the density and then vary a through a Feshbach resonance (Bloch et al. 2008). 
Then we can move from the BCS regime (kpa — > — ), where A oc exp(—l/kp\a\) and 
fx oc e F (Leggett 2006) , through unitarity at the threshold for two-body binding, where 
| a, | 1 — > 0, and into the deep BEC regime (kpa — > + ), where /i oc —h 2 /2ma 2 . The 
latter shows that the energy per particle is half the binding energy of the molecules 
in the deep BEC limit as it should be. Also, one can show that in the BEC limit, 
A oc ep /(na 3 ) 1 / 3 . For the pair wave function of (271, this implies that £' oc a. The 
bound state wave function in figure[3]is then exactly of the expected form, e~ r ^ a /r, and 
becomes a delta function for very tightly bound pairs (a — > 0). Importantly, keeping 
a constant and varying n is another way that one can explore the crossover. This 
can be seen by realizing that in the simple Leggett model, the dimensionless control 
parameter is kpa oc an 1 / 3 . This observation will be important for the application 
discussed in Section 

The crossover is then the sequence of structures occurring as the interaction is 
varied from weak attraction describable in BCS-theory, via stronger attraction barely 
binding pairs, to strong attraction with substantial binding energy of pairs. When 
pairs can bind it is energetically advantageous to place all bound pairs in the same 
lowest mean-field level for pairs, thus forming a BEC of these bosons made of two 
fermions that exists only in the confinement of the external field. The transitions in 
this sequence is driven by adiabatic or sudden changes of the interaction or equivalently 
of a, resulting in corresponding dynamic evolution of the structure. This can either 
be evolution towards an equilibrium or by establishing quasi-stationary coherent 
oscillations between the initial state of paired atoms to the structure of a molecular 
BEC (Pethick & Smith 2008, Bloch et al. 2008, Giorgini et al. 2008). 

The unitary region is of considerable interest since \a\ — »• oo means that the system 
effectively loses a scale, and the divergent scattering length can thus not be used as 
expansion parameter when calculating the energy of the system (Bertsch 1998, Giorgini 
et al. 2008). In this situation, the structure has to be independent of a, and hence 
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is referred to as the universal limit. For a Fermi system, we thus expect the energy 
to be proportional to the only scale in the problem: E/N = |(1 + /3) 1 m F , where (3 
(sometimes given as r\ = 1 + j3) is a parameter that must be determined numerically. 
Monte Carlo approaches (Carlson et al. 2003, Astrakharchik et al. 2004, Chang 
et al. 2004) give a value of f3 = — 0.58±0.01, which indicates considerable more binding 
than the BCS mean-field result of (3 = —0.41. Analogous universal relations can be 
defined also in the imbalanced case (Giorgini et al. 2008) where the two internal states 
of the fermions are unequally occupied. The bottom line is that any Fermi system is 
expected to exhibit similar universal behavior provided a is large compared to other 
length scales. 

In atomic physics, the simple homogeneous version of the gap equation in Eq. (22 1 
is most commonly used. The assumptions are low density, negligible effects from the 
trap, and consequently plane wave states. In this dilute limit distances are large, 
momenta small, and only binary encounters are expected to contribute to the physics. 
These conditions are met in a typical atomic physics setup, where the density, n, is 
fixed and the scattering length, a, is varied. 

Low density and plane waves are usually not applicable in finite nuclear systems, 
where the basic interaction in principle is fixed by the nucleon-nucleon potential. The 
bare (physical) neutron-neutron scattering length (<~ —18 fm) is an order of magnitude 
larger than the range of the nuclear force (~ 1.4 fm), see Tab. ([I]). Thus, there is 
hope that a system of low density neutrons would show universal behavior and be 
described by something resembling a BCS-BEC crossover model. The approximately 
infinite neutron matter that one might expect to find in a neutron star therefore seems 
like an obvious candidate for applying crossover theory to nuclear systems (Lattimer 
& Prakash 2004). 

However, the need for effective nucleon-nucleon potentials complicates the 
problem. This leads us back to the problem of separating mean-field and residual 
interactions for nuclei. To illustrate, we note that if one uses BCS theory in infinite 
nuclear matter with a realistic bare neutron-neutron potential, then one finds a value 
that is much smaller than both the measured gap and the value from the otherwise 
successful semi-empirical mass formula (Fetter & Walecka 1971). Of course, finite 
systems yield different results, but for large nuclei one should expect fair agreement 
with this limit. An effective neutron-neutron interaction can be used to give a better 
gap, but this is then only applicable in a stricly defined limited context in a local 
area of the nuclear chart and within a particular approximation scheme (such as BCS 
theory) . 

The unitary Fermi gas and the BCS-BEC crossover dynamics that can be studied 
with cold atomic Fermi gases also plays an important role in elucidating the nature 
of the so-called 'pseudogap' regime which has been proposed in connection with high- 
temperature superconductivity (Chen et al. 20056, Chen et al. 2009). Theoretical 
studies find a pseudogap phase above the critical temperature for superfluidity due 
to preformed pairs (Randeria 1998, Janko et al. 1997, Yanase & Yamada 1999, Perali 
et al. 2002, Bruun & Baym 2006, Barnea 2008). However, these preformed pairs 
are not condensed due to large fluctuations and the underlying Fermi statistics 
of the particles. These predictions have been confirmed by ab-initio Monte Carlo 
calculations (Magierski et al. 2009, Magierski et al. 2011) and can also be obtained 
using quantum cluster expansions (Hu et al. 2010). A number of recent experiments 
(Stewart et al. 2008, Gaebler et al. 2010, Perali et al. 2011) provide strong support for 
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preformed pairs and for the pseudogap theory (Randeria 2010). 
4. The Atom Gas from a Nuclear Perspective 

In order to meaningfully discuss the similarities and differences between the atomic 
and nuclear physics systems it is imperative to understand the relevant quantum 
states involved in the two situations. Therefore we will now discuss the much studied 
two-component Fermi gas from atomic gas physics. Bosonic states in nuclei will be 
addressed later in relation to a-particles. Then we will compare and contrast the 
interaction in the atomic physics systems with that typically employed in nuclear 
physics in the framework of the (nuclear) shell model. We shall continue to emphasize 
corresponding degrees of freedom in nuclei and cold atomic gases. Also for convenience 
we have collected in Table ([I]) a number of typical values for characteristic key 
quantities in the two systems. 

4-1- Two-component Atomic Fermi Gases 

By now the production of degenerate Fermi gases has been achieved by many 
experimental groups around the world (Bloch et al. 2008, Giorgini et al. 2008). 
However, this was initially a considerable challenge since the evaporative cooling of a 
single species of fermions is hindered by the Pauli exclusion principle since they are 
basically non-interacting. This means that the gas cannot equilibrate and lower its 
temperature by expelling the fastest moving atoms. If one instead prepares a gas with 
a population of two different species of fermions, then low-energy s-wave interactions 
are no longer forbidden and degeneracy can be reached (Ketterle & Zwierlein 2008). 
Typically the Fermi temperature, Tp — h 2 k F /2mks, is around Tp ~ 1/iK, whereas 
experiments can cool the gas to T < 200nK, thus T/Tp < 0.2 so the degenerate regime 
is truly reached. The typical density of atoms in the trap is around 10 13 cm~ 3 . This 
translates to an average interatomic distance of about 500 nm, which is much bigger 
than atomic sizes (of order 0.1 nm) and inter-atomic potential ranges of order 1-10 nm. 
In contrast, nucleons in nuclei have average interparticle distances of the order ~ 1 fm 
which is the same order as the nucleon size (see Table ([T])) (Bohr & Mottelson 1975), 
and nuclei are never dilute in the sense of atoms in an ultracold trapped gas. The 
criterion for diluteness in the atomic case is usually no? <C 1 (away from Feshbach 
resonances). In nuclei, na? is of order one or larger. This is one reason to be careful 
in comparing the two systems. 

The atoms must be confined by an external field, a trap, and two distinctly 
different degrees of freedom are present, i.e. internal ones related to the atom itself 
and external ones related to the atoms in the trap. Correspondingly, two different 
energy scales emerge, i.e. the spacing of (external) trap eigenstates and the energy 
difference of the (internal) atomic states. We will now discuss these in detail for the 
relevant experimental setup. 

The internal atomic quantum numbers are those of the standard hyperfine 
interaction, obtained from the coupling of the angular momentum (orbital, Z, and 
spin, S) of the electrons and the nuclear spin, 7, into total spin F and projection 
mp. The complete set of quantum numbers is (J , J, F,mp). In a magnetic field 
only mp is conserved. In the limits of high and low field strengths there are good 
asymptotic quantum numbers (Pethick & Smith 2008). Systems with two components 
corresponding to population of the lowest two hyperfine states (of two different given 
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projections mp) have been extensively studied (Ketterle & Zwierlein 2008), and 
recently also three-component systems (three different mp values) have been realized 
in experiments (Ottenstein et al. 2008, Huckans et al. 2009, Nakajima et al. 2010). 
The hyperfine energy splitting for different mp due to the magnetic field is of order 
10~ 6 eV for the alkali atoms typically used in experiment (Pethick & Smith 2008). 
This is much smaller than any electronic transition in the atoms. 

The external quantum numbers depend on the trap which often is harmonic 
and we shall use those of the harmonic oscillator, which for the isotropic case is 
(nlm). Typically, the oscillator (trap) length b = \Jhj Am^u (where mjv is the 
nucleon mass, ui the trap angular frequency, and A the mass number of the trapped 
species) is of the order of in a large trap. This can be translated into an energy 
huj = 4.18- 10 _11 A _1 (l/im/6) 2 eV. We immediately see that the trap level spacing hoj 
is much smaller than the hyperfine splitting of order /zeV. In total, the two-component 
Fermi gas in a trap therefore has two different internal states, mp and m' F , and each 
atomic state is in addition to the m^-value associated with external quantum numbers, 
{nlm). 

If we assume a distribution of identical fermions occupying the (lowest two) 
hyperfine states, the lowest oscillator levels would be occupied for each mp depending 
somewhat on the temperature. Inelastic collisions can in principle change the hyperfine 
projections of a pair of atoms. These processes are, however, suppressed and the two- 
component Fermi gas is therefore in practice stable with a fixed number of atoms 
in each hyperfine state (Ketterle & Zwierlein 2008). The hyperfine states can be 
considered as frozen degrees of freedom when we describe the system. The internal 
and external degrees of freedom are uncoupled, and no spontaneous process exist that 
changes the internal quantum numbers on the timescale of up to seconds of a typical 
experiment (Ketterle & Zwierlein 2008). This decoupling is what allows one to treat 
separately these two sets of degrees of freedom and use product wave functions. 

Beside the one-body external field also two-body interactions between the atoms 
are present, thereby introducing another energy and length scale into the problem. The 
central part of the two-body interactions between the atoms is of the van der Waals 
type with characteristic length scale of 10-100 Bohr radii (0.5-5 nm) and typical low- 
energy s-wave scattering lengths can be 100-1000 Bohr radii and of both signs (ignoring 
for now the interesting case of Fcshbach resonances) . The average distance (500 nm) 
in the cloud estimated above is larger than these interaction lengths. This would 
imply that two-body interactions are weak. The ratio of the interaction energy, Ej to 
the energy of the trap, Eho can be estimated for TV atoms in a harmonic oscillator 
mean-field potential of length b and a ^-interaction of strength given by the scattering 
length a, i.e. Ei/E HO ~ iV 1/6 a/6 (Pethick & Smith 2008). If we take typical numbers, 
a ~ 10 — lOOnm and 6^1 — 10/im, we therefore see that with a/b ~ 0.1 — 0.001 the 
ratio will be small even for millions of particles. The effect of interactions on the 
ground state away from resonance is therefore usually small. 

The two-body interaction may, however, be varied experimentally by utilizing 
Fcshbach resonances resulting in anything from attraction to repulsion corresponding 
to scattering lengths in the interval [— oo;oo] (Chin et al. 2010). The Feshbach 
resonance technique exploits two internal atomic two-body states of different magnetic 
moments, which means that the relative energy of the two states can be controlled by 
an externally applied magnetic field. We illustrate the situation in figure [4] When the 
scattering energy in the in-coming (open) channel is resonant with a bound state in a 
closed channel of different magnetic moment, the scattering length blows up and passes 
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Figure 4. Illustrative plot of the physics of a Feshbach resonance. The lower 
(blue) potential curve indicates the open (in-coming) channel of two atoms in a 
given spin state, while the upper (black) curve is the closed channel that is in a 
different spin state from the open channel. The dotted horizontal line indicates 
the relative energy of the incoming atoms. The horizontal (red) line in the closed 
channel indicates the presence of a bound state. By tuning an externally applied 
magnetic field the open and closed channel potential curves can be moved relative 
to each other. In A) the bound state is above the relative energy, while it is goes 
below in B). The characteristic divergence of the scattering length, a, happens 
when the bound state energy coincides with the reletive energy of the incoming 
atoms, and a changes sign when going from situation A) to B). 

through infinity as the magnetic field is varied. One thus has laboratory control of 
the effective interaction of the atoms in the (open) scattering channel. This situation 
is connected with the number of bound states of the two-body interaction. The point 
at which the scattering length diverges is where another bound state appears (or 
disappears) in the potential. Note that Feshbach resonances which are properties of 
the two-body potentials occur for all atomic systems, irrespective of whether they 
are bosonic or fermionic species. However, the s-wave interaction can of course 
only occur between two bosons, a boson and a fermion, or two different internal 
states of fermions. For a single-component Fermi gas (fully polarized) the dominant 
interactions are usually in the p-wave channel. We also note here that there is a slight 
abuse of common language in scattering theory where the word 'resonance' usually 
denotes a feature (typically a peak) in the scattering cross section. In the realm of 
ultracold atomic gases, a Feshbach 'resonance' most often refers to the divergence in 
the scattering length as a certain external (usally magnetic) field is varied. The cause 
of this behavior of the scattering length is of course associated with a resonance in the 
traditional sense between an open and a closed channel as explained above. 

A common parametrisation of the s-wave scattering length variation as the 
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magnetic field is tuned around a Feshbach resonance takes the form 

a(B) = a bg (l - ^^-) , (28) 

where Bq is the magnetic field position corresponding to a resonance at zero energy, 
AB is called the resonance width, and ab g is the background value of the scattering 
length away from the resonance. AB and a^g can be of both signs in general. Values 
and parameters for Feshbach resonances in different atomic systems can be found 
in (Chin et al. 2010). Correlations may increase tremendously when the scattering 
length diverges as the resonance is approached. Different features like BCS-BEC 
crossover and universality as in Efimov physics (Section [6| then arise with a demand 
to describe the correlations. This can be compared with the nuclear physics case where 
the predetermined strong nuclear interactions causes a highly correlated ground state 
configuration. Mean-field effects are still prominent but correlations often contribute 
significantly. 

We now estimate the ratio of interaction energy to mean-field energy in a nucleus 
in order to compare to the atomic physics ratio found above. We use a crude 
harmonic oscillator value for the mean-field part and a phenomenological estimate 
for the residual pairing interaction. The mean-field energy is Nep ~ NhujjGN) 1 / 3 , 
where 6f is the Fermi energy, and we have Hoj = 41/\ZiV MeV (see Table. (111)). The 
energy gain from the pairing interaction is proportional to the single-particle level 
density gp « N/10 MeV -1 at the Fermi energy times the square of the pairing gap 
A w 12 /VN MeV (Bohr & Mottelson 1975). Thus the ratio is about 1/(5V 5 / 6 ) which 
is small and further decreases with the number of nuclcons. In terms of this estimate, 
the nuclear and atomic system are only comparable in the weak-coupling limit where 
\a\ — > 0. When the atomic scattering lengths are large, the nuclear and atomic systems 
are basically incomparable, i.e. the fundamental interactions are completely different 
implying that the physics also is different. We note that the weak-coupling limit is not 
very interesting in nuclei since there is a critical interaction strength for superfluidity 
in finite Fermi systems. However, this does not imply that the methods and techniques 
cannot be carried from nuclei to atomic gases. Advanced nuclear models can deal with 
correlations in both fields. 

In studies of trapped (two-component) atomic Fermi gases one can often replace 
the discrete one-body oscillator levels by the homogeneous spectrum, essentially 
switching to continuous momentum space. This is because the oscillator level spacing 
is very small and the spectrum can be well approximated as a continuum. This also 
holds for strong interactions, where local density or Thomas-Fermi approximations 
often are applied (Pethick & Smith 2008, Pitaevskii & Stringari 2003, Giorgini 
et al. 2008). 

The temperature scale of ultracold atomic gas experiments is in the range T ~ 
10 — 100 nK which is extremely small but nevertheless non-zero. The corresponding 
energy ksT ~ 10~ 12 — 10~ n eV is larger than the typical atomic trap spacing. Thus, 
the spectrum is approximately continuous and sums can be converted to appropriately 
weighted integrals in thermodynamic considerations. If the trap becomes tighter (small 
b) this approximation is no longer expected to be accurate as the level spacing of the 
trap states goes as Huj cx b~ 2 . Then one must explicitly use the discrete spectrum 
and the quantum numbers of the harmonic trap. In this respect, tight traps with 
very few particles is the case most closely related to nuclei (Serwane et al. 2011, Ziirn 
et al. 2012). Alternatively, a very deep periodic optical lattice potential can have 
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single sites that are approximately harmonic and contain only a few particles (Bloch 
et al. 2008), again closer to the situation in small nuclei. The notion of temperature 
is, however, somewhat ill-defined for isolated systems like self-bound nuclei where 
the energy is conserved. The temperature is only useful in connection with average 
quantities in sufficiently excited nuclei where the level spacing is relatively small (Bohr 
& Mottelson 1975, Siemens & Jensen 1987). 

The tunability of atomic interactions through Feshbach resonances is an 
important feature of experiments that enable studies of physical system over broad 
regions of parameter space. In particular, one is forced to think more in terms of 
generally applicable theories. In nuclear physics a number of scales are intrinsic to 
the problem and do not change, even when considering stellar environments such 
as neutron stars. Typical values are listed in Tab. (JlJ. A prominent example is the 
neutron- neutron scattering length which is a ~ —18 fm. This is an order of magnitude 
larger than the typical interaction range in nuclei of ro ~ 1.4 fm. This would suggest 
that neutron matter is in a universal regime, i.e. that the physics does not depend 
on the specific details of the inter-neutron potential only on the scattering length 
and the density of the neutron matter. The dimensionless measure of the strength of 
the interaction in this situation is kpa with kp oc n 1 / 3 the Fermi wave vector for a 
homogeneous system. This is similar to a two-component atomic Fermi gas close to a 
Feshbach resonance where a diverges. However, the correction to the scattering length 
description which is given by the effective range is intrinsically large in nuclei, r e ~ 2.7 
fm (Siemens & Jensen 1987). So when kpa is large so is kpr e . A theory that takes 
the effective range correction into account was discussed in (Schwenk & Pethick 2005) 
while a recent numerical study can be found in (Gezerlis & Carlson 2010). Problems 
with non-universal corrections also arise in connection with bosonic few-body systems, 
when the range of the potential, ro, is neglected and only the scattering length is kept, 
i.e. ro/a or r e /a are small, as is often done in effective field theory calculations of few- 
body properties (Deltuva k. Lazauskas 2010). 

In the cold atomic gases there is, however, a handle also on the relation between 
a and r e near a Feshbach resonance. This stems from the fact that the effective range 
can be related to the width of the Feshbach resonance, AB (Chin et al. 2010). In fact, 
in zero-range models, r e (B) on resonance (|o| = oo) is given by 

2h 2 

reO = XWa~> (29) 

with m the atomic mass and A/x the difference in magnetic moment between the 
coupled channels that cause the resonance (Petrov 2004, Jonsell 2004). The effective 
range can thus be reduced by working with broad resonances. Many cold atomic Fermi 
gas experiments have aimed for the universal physics as a diverges and have worked 
with broad resonances where background scales are small (Bloch et al. 2008, Ketterle 
& Zwicrlein 2008). Also the wide resonances are easier to work with from an 
experimental point of view as they require less resolution on the applied magnetic 
field. However, very recent experiments have now started utilizing narrow resonances 
(Gross et al. 2010). 

An aspect of recent interest for experiments is the ability to tune the scattering 
length to zero at 5 = Bo + AB in the parametrisation of Eq. (28) (Roati 
et al. 2007, Pollack, Dries, Junker, Chen, Corcovilos & Hulet 2009). This was 
used to reduce interactions in atom interferometry (Fattori, DErrico, Roati, Zaccanti, 
Jona-Lasinio, Modugno, Inguscio, & Modugno 2008) and to probe magnetic dipolar 
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interactions (Fattori, Roati, Deissler, DErrico, Zaccanti, Jona-Lasinio, Santos, 
Inguscio & Modugno 2008). In zero-range models of Feshbach resonances, the effective 
range is a function of the magnetic field, r e (B), and in fact diverges when a — > 
(Zinner 2009), a fact that also holds for a model potential like the attractive 1/r 6 
potential which is the outer part of the atomic van der Waals interaction (Gao 1998). 
The only remaining scales are then ab g and r e o from Eq. (29). For a bosonic system 
this can lead to strong modification of the stability properties around the resonance 
(Zinner & Th0gersen 2009, Th0gersen, Zinner & Jensen 2009, Zinner 2011). For atomic 
Fermi gases the regime where a — > remains largely unexplored which is due to the 
fact that the experimental use of broad resonances dominates the field and makes this 
non-interacting limit less interesting at the moment. As experiments become more 
sofisticated it is, however, clear that one should include the non-universal corrections 
from a non-zero effective range in the description of experimental findings. 



4-2. Interaction in a Shell Model Picture 

We now consider the interaction in the atomic gases and compare the language to the 
treatment of interactions in the standard nuclear shell model through two-body matrix 
elements. This will also lead us to a discussion of how to map the two-component Fermi 
gas with and without trap into the nuclear system of a single nucleon species. We also 
comment on the possibility of using the isospin formalism to describe multi-species 
Fermi systems. 

As mentioned, the interactions of typical two-body terms in atomic gases and in 
nuclei are vastly different (see Tab. ([!])). The atomic interactions are of ranges that 
are several orders of magnitude smaller than the average distance in dilute gases. For 
nuclei it is quite the opposite. Here the range is a few Fermis which is similar to the 
radius of the nucleus. The observation of pairing gaps and superfluidity in atomic gases 
demonstrates that interactions are very important, so in this respect atomic gases are 
similar to nuclei. However, we will now address a very fundamental difference that 
can cause confusion. 

In the atomic gases, the two-body interaction originates from the van der Waals 
force and is a real physical interaction (Weiner et al. 1999, Kohler et al. 2006). The 
Hamiltonian is divided into an external one-body piece given by the trap and the 
two-body interaction from the atom-atom collisions governed by the van der Waals 
interaction. In the nuclear case, this kind of splitting is not provided a priori as 
the nucleon-nucleon interaction is all there is. However, in practice most methods 
find it convenient to split the interaction into a mean-field and a residual two-body 
interaction. The mean-field is chosen to reproduce average bulk properties such as 
energies and shell closures, whereas the residual interaction should reproduce features 
beyond mean-field such as collective states, supcrfluid properties, etc. 

In practice, this splitting for nuclei is rather arbitrary and context 
dependent, because the basic as well as the effective nucleon-nucleon interaction is 
phenomenologically adjusted to selected observables. With a chosen interaction the 
split is seemingly optimized with the self-consistent mean-field (Ring & Schuck 1980). 
However, the very choice of effective interaction depends on both the selected 
observables and the chosen method of computation e.g. resulting in a limited Hilbert 
space (Siemens & Jensen 1987). The residual interaction is then whatever is left 
over from any given mean-field computation. These are essential problems of nuclear 
physics that have plagued the field. 
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Table 1. Estimates of typical parameters and ratios of interest in nuclear and 
cold atomic gas systems in natural units of MeV and fm for nuclei, and eV and nm 
for atoms. The bosons in the table under nuclear physics are effective constituents. 
For cold atomic gases we merely cite a few alkali species of recent interest, n is the 
particle density, rrj is the characteristic distance of the interparticle force, which is 
the van der Waals length for cold atomic gases. The scattering length is given for 
both neutron-neutron (nn) and neutron-proton (np) for nuclei. N is the number 
of atoms (with two spin states for the two-component gas considered in the table) 
and A is the number of nucleons. b is the harmonic trap size for atoms and the 
system size for nuclei, ftw is the trap level spacing of the outer confinement for 
atomic gases, whereas for nuclei it is the level spacing in a mean-field oscillator 
model (Fetter & Walecka 1971). E HO ~ NhuiiQN) 1 / 3 is the total mean-field 
energy for atoms and Ejjq ~ Nhui(3A/2) 1 ^ 3 for nuclei, ejj- = fr 2 /mr 2 , is the 
characteristic kinetic energy of the interparticle potential and ej n t is the energy 
difference between the internal degrees of freedom. For atoms this is the hyperfine 
scale AEhf, whereas for a single nucleon system with equal spin up and down 
components it is rouhgly the mean-field splitting hu>. Ej is the average value of a 
typical residual interaction for nuclei and for atoms we use Ej/Ejjq ~ A rl ' 6 a/6 
as in the text. A is the typical BCS pairing gap. The numbers for atoms given 
in the latter part of the table are for N = 10 , expect for the Fermi energy, Ep, 
which is given by Ep = (37r 2 ) 2 / 3 7i 2 n 2 / 3 /2m where m has been taken as the 6 Li 
mass for atoms and the nucleon mass for nuclei. 



Quantity 


Nuclear Physics 


Cold Atomic Gases 


Fermions 


neutron, proton 


6 Li, 40 K, etc. 
7 Li, 85 Rb, etc. 


Bosons 


dcutcron (np), a 


n 


0.10-0.15 fm- 3 


10 10 - 10 13 cm~ 3 


d = n~ 1 / 3 


~ 2 f m 


500-5000 nm 


ro 


~ 1.4 fm 


0.5-5 nm 


a (non-resonant) 


-18 fm (nn) 5 fm (np) 


(±)10-100 nm 


n|a| 3 


10 1 - 10 3 


10~ 6 - 10~ 2 


A/N 


10-100 


10 3 - 10 6 


b 


l^AVSfm 


1-10 (an 




4,1 /VA MeV 


< 10- 11 eV 


Eho 


47A 5 / 6 MeV 


< 10" 3 eV 




10 MeV 


10- 7 - 10- 5 eV 


tint 


~ hw 


AE hf ~ 10" 6 eV 


Ei 


5-15 MeV 


< 10~ 4 eV 


Ep 


~ 43 - 56 MeV 


HT 12 - HT 10 eV 


A 


0.5-2 MeV 


< 0.5t F 


E z /Ne K 


0.02-0.2 


0.01-1 


Ej/E H o 


0.006-0.04 


0.1-0.001 


A/ep 


< 0.05 


< 0.5 



In the following we will point out some of the problems that different choices 
of mean-field and residual interactions can cause. In the atomic case the choice is 
in this sense much more pure. At most it can be debated where to put mean-field 
contributions from the two-body interaction. Great care must be exercised to provide 
a useful and meaningful comparison which can be used to transfer techniques and 
concepts. 

Let us now discuss the atomic interactions from a nuclear physics perspective. 
First we argue that s-waves of the relative motion of two scattering atoms are the 
dominant two-body term for the two-component gas. To distinguish from the external 
trap quantum number I we name the two-body relative orbital angular momentum l r . 
Recall that the gas is extremely dilute with average atom-atom distance much larger 
than the interaction range. The two-body wave function approaches zero as r 2lr , 



CONTENTS 



29 



where r is the relative distance. For l r ^ there is a very small probability of finding 
the particles at small r-values, whereas for l r — the probability is non-zero. For very 
short-range interactions, it is therefore a good approximation to restrict interactions 
to the s-wave channel. 

The l r — two-body states are symmetric, and the Pauli principle can only 
be satisfied by using an antisymmetric combination of the other parts of the wave 
function. If the system is fully polarized with all atoms in a single hyperfine state, then 
s-waves would be Pauli forbidden. The atoms therefore essentially do not interact, 
only very weakly through odd parity Pauli allowed relative states, i.e. dominantly by 
p-waves. On the other hand, two atoms in different hyperfine states can interact much 
stronger through the allowed relative s-wave interaction. 

In nuclei the situation is quite different as the nucleon-nucleon interaction has a 
hard-core part which generates terms of all multipolarities where the strength diminish 
slowly for increasing l r (Bohr & Mottelson 1975, Siemens & Jensen 1987). Furthermore 
the spin-dependence is quite complicated and spin-dependent terms are strong in 
nuclei. This is why a full analogue with dilute gases of atoms is only appropriate 
for low density nuclear matter (ideally low-density neutron matter to avoid Coulomb 
interaction) or perhaps in very neutron-rich or halo nuclei. We will discuss some 
examples in the following sections. 

4-3. Mapping atomic and nuclear degrees of freedom 

The shell model has been extremely successful in nuclear (structure) physics. First 
a mean-field is generated, self-consistently and often phenomenologically adjusted, 
as a basis for a more detailed description. The preferred non-sclf-consistent 
phenomenological choice is the so-called Woods-Saxon potential (Bohr & Mottelson 
1975), which provides bound state single-particle properties similar to that of a 
harmonic oscillator. For nucleons in the mean-field, a sizable spin-orbit coupling 
is present with opposite sign compared to the spin-orbit term for electrons in atoms. 
Due to the strong spin-orbit coupling in nuclear physics, it is customary to work in the 
coupled representation where orbital angular momentum I and nucleon spin s = 1/2 
is coupled to j = I ± 1/2 (Fetter & Walecka 1971). 

For atomic physics applications of nuclear methods this coupled basis is not the 
most convenient choice. To spell out the precise analogies for a two-component atomic 
system we relate the external atomic trap degrees of freedom with the nuclear orbital 
motion, and the internal atomic hyperfine states with the two spin projections of the 
nucleon. Since the external and internal atomic degrees of freedom are effectively 
decoupled, an uncoupled product basis of nuclear orbital and spin motion is the 
closest analogy. This essentially corresponds to a harmonic oscillator without spin- 
orbit coupling. 

To be specific we consider a two-body matrix element in the interacting nuclear 
shell-model with the interaction relevant for the two-component atomic gas. First we 
have to specify the mapping of atomic states to nuclear states. For this we choose 
to work with a single species of a spin 1/2 nucleon. The obvious thing is to map 
the two hyperfine states of the two-component atomic gas to the two spin projections 
of the nucleon. The nuclear two-body spin wave functions can be either singlet or 
triplet (total spin of or 1) which corresponds to atomic antisymmetric or symmetric 
(internal) hyperfine combinations, respectively. 

The atomic two-body interaction acts in the s-wave (external) channel on opposite 
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hyperfine (internal) states only. This maps to s-wavc nuclear interactions between 
opposite nucleon spin projections only. Therefore, the s-wave atomic interaction in 
nuclear physics notation should be expressed as V r atom (r) = Ps=oV2{r)Ps=o, where 
V^(r) is the spatial part of the interaction in terms of the relative coordinate r 
and Ps=o = (1 — o\ ■ <?2)/4 is the projection onto the spin singlet (antisymmetric) 
component of the two-body state. One can now use the standard methods of angular 
momentum coupling to transform between coupled and uncoupled representations and 
thus calculate any matrix element needed. An explicit calculation and transformation 
for the zero-range interaction can be found in (Ozen & Zinner 2009). 

An important point in the mapping above is that the two hyperfine atomic states 
in a two-component Fermi gas are mapped to a single species of nucleon with each 
nucleon spin direction mapped to one of the two hyperfine states. Of course there are 
both protons and neutrons in nuclei and one could ask whether this has interesting 
analogies in atomic physics. The first obvious system to look for is that of trapped 
gases with two species of fermionic atoms, each with two internal hyperfine states. 
In this case the s-wave interactions between identical atoms in different hyperfine 
states would remain but between unlike atoms all combinations of hyperfine states 
would have finite interactions. This could resemble the situation in nuclear physics 
where neutron-neutron, neutron-proton, and proton-proton s-wave interactions are all 
allowed including extensively studied pairing interactions of different types (Dean & 
Hjorth- Jensen 2003). However, the neutron and proton interactions should then not 
be related as in nuclear physics but translated into independent atomic interactions 
each connected to their own individual hyperfine states. Recently such systems have 
been seriously explored in alkaline-earth gases with the purpose of use for quantum 
computation (Daley et al. 2008, Gorshkov et al. 2009), and for testing S£/(iV)-niodc-ls 
for large N (Gorshkov et al. 2010, Xu 2010). Experiments are also on-going (Fukuhara 
et al. 2009, De et al. 2009). Nuclear physicists often use isospin as a quantum number 
representing the similarity of neutrons and protons in strong interactions (Siemens & 
Jensen 1987). Isospin is a measure of the symmetry of the wave function resulting from 
the permutation symmetry of neutrons and protons in the Hamiltonian. Therefore, it 
is not obvious that isospin has an analogue in the atomic systems. This is because the 
scattering properties of atoms depend sensitively on species and hyperfine states and 
would therefore most likely strongly break this symmetry. The alkaline-earth gases 
above are, however, an example where the spin of the atomic nucleus can be mapped 
onto isospin (Gorshkov et al. 2010). 

Another possibility is a two-component system of bosons where particles in the 
same hyperfine state are allowed. The three interactions between pairs in the same and 
in different hyperfine states should then be similar to an isospin 1 system. Also a three- 
component fermionic atomic system has three possible s-wave interactions (Ottenstein 
et al. 2008, Wenz et al. 2009) and could resemble an isospin 1 system. However, it must 
be kept in mind that generally the scattering lengths arc different between different 
hyperfine components. Treating them as equal can be a good approximation at times 
and may help understand the symmetries of the system. Of course, a full quantitative 
analogy requires that one takes the different scattering lengths into account directly. 

A concrete example of a fruitful interaction between nuclear physics methods and 
the physics of cold atomic gases is the pairing gap (Bohr & Mottelson 1975). This can 
be defined independently of any particular theory, as 



A(N) = 



E(N + 1) - 2E{N) + E(N - 1) 
2 



(30) 
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with E(N) the ground-state energy of the -/V-body system. Extrapolation of the 
pairing gap in nuclear physics from binding energy differences as in ( 30 ) presupposes 
that the underlying structure is unchanged by adding or subtracting one nucleon. This 
assumption is not always well fulfilled as for example deformation may change rapidly 
when closed shell configurations are approached. A way to remove this problem to 
second order in such smoothly varying effects (such as deformation) is to use instead 
the double difference (Jensen & Hansen 1984) 

- [E(N - 2) - 3E{N - 1) + 3E(N) - E(N + 1)] (31) 

This is usually an improvement although a slightly larger range of nuclei are employed. 
For cold atomic gases such a procedure is typically not necessary as the expression in 
(30 1 is sufficiently accurate. 

The gap is a measurable quantity, both for nuclei and for ultracold Fermi 
systems. In the nuclear setup, there are many advanced many-body methods that 
go considerably beyond BCS theory to calculate the gap. As ultracold gases become 
smaller and finite size and small particle numbers begin to be important, these 
methods can be directly transfered from nuclear physics by simple re-ordering and 
re-interpretation of models and interactions as we have discussed above. 



4-4- Fermionic few-body systems 

A particular promising venue where nuclear physics and ultracold atoms can have 
large overlap is in the study of fermionic systems with a limited number of particles. 
This is a natural condition for the physics of light nuclei with intrinsically fermionic 
nucleon constituent. Usually the number of atoms in a cold atomic gas experiment 
has been of order 10 3 — 10 6 and a theoretical approach based on few-body methods is 
obviously doomed. However, recently microtraps have been built that are capable of 
trapping very few atoms (1 — 10) at very low temperatures (Serwane et al. 2011, Ziirn 
et al. 2012). Furthermore, it is now possible to address single lattice sites in an optical 
lattice and produce and probe only a few atoms on such a site (Bakr et al. 2010, Sherson 
et al. 2010, Weitenberg et al. 2011). These developments provide hope that aspects of 
nuclear physics can be simulated with cold atomic gas setups in a very direct manner. 

In the meantime, a question of immense theoretical interest for fermionic two- 
component systems was posed back in 1999 by George Bertsch and concerns the 
nature of the ground-state of a system of fermions with two internal states when the 
interaction between fermions in different internal states have a large scattering length 
( going to infinity, i.e. the unitarity limit). The relevance of this question for atomic 
gases has been outlined above, whereas for nuclear physics it is relevant for instance 
in neutron matter where the scattering length in the singlet channel (opposite spin 
direction) is large compared to the range of the nucleon-nucleon interaction. This 
question has sparked immense activity both for large and for small systems (a recent 
review can be found in (Giorgini et al. 2008)). Here we are interested in the connection 
to nuclei and to cold atom experiments with a limited number of particles. We will 
therefore focus almost exclusively on small particle number. 

Early studies of this interesting question in the homogeneous case have already 
been discussed in Section |3.3| and here we focus on the case where the particles are 
confined by an external harmonic potential. An analogy to nuclei can be made by 
considering a mean-field potential where nucleons are confined to interact through the 
residual interaction. 
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Table 2. Energies (in units of hw) calculated with the SMMC method for a 
trapped fermion gas with scattering lengths a/b = 11 (BEC) and a/b = —1.0 
(BCS) for different particle numbers N. The statistical uncertainty is given 
in parenthesis. HOSD denotes the non-interacting energies. From (Zinner 
et al. 2009). 



N 


HOSD 


BEC 


BCS 


N 


HOSD 


BEC 


BCS 


2 


3 


1.72(3) 


2.49(3) 


12 


32 


20.7(2) 


27.23(2) 


3 


5.5 


3.9(2) 


4.84(2) 


13 


35.5 


24.0(2) 


30.21(2) 


i 


8 


4.96(3) 


6.84(4) 


14 


39 


26.0(1) 


33.16(2) 


5 


10.5 


7.1(3) 


9.14(2) 


15 


42.5 


29.1(1) 


36.08(5) 


6 


13 


8.11(7) 


11.08(5) 


16 


46 


30.9(1) 


39.03(2) 


7 


15.5 


10.6(2) 


13.28(5) 


17 


49.5 


33.6(2) 


41.98(2) 


8 


18 


11.58(5) 


15.21(4) 


18 


53 


35.7(1) 


44.89(2) 


9 


21.5 


14.8(2) 


18.30(4) 


19 


56.5 


39.0(1) 


47.85(2) 


10 


25 


16.34(6) 


21.27(3) 


20 


60 


40.8(1) 


50.75(2) 


11 


28.5 


19.3(2) 


24.28(2) 


21 


65.4 


44.5(2) 


54.50(2) 



The early seminal study of Busch et al. (Busch et al. 1997) demonstrates that 
for two trapped fermions in opposite internal states the ground-state energy is 2huj 
(see (Zinner 2012) a discussion of this model in both two- and three-dimensional 
traps and for a review of the relevant theoretical work and experimental support). 
An important benchmark for the two-component fermionic few-body problem in a 
harmonic trap is the exact solution for three particles at unitarity where the ground- 
state energy is 4.27/iw (Werner & Castin 2006a, Werner & Castin 20066). These 
results was followed by several numerical studies of three and four trapped fermions 
(Stetcu et al. 2007, Alhassid et al. 2008, von Stecher & Greene 2007), and for up 
to 30 fermions (Chang & Bertsch 2007, Blumc et al. 2007). All these studies found 
agreement with the exact result for three particles and furthermore established that 
five fermions have an energy of 5fiw to within a few percent at unitarity. We note that 
these nice universal results at unitarity break down if the two kinds of fermions have 
different masses at which point finite-range corrections must be taken into account 
(Blume & Daily 2010). 

A number of studies have pursued the energy and structural properties of 
fermionic systems as the scattering length is varied across the Feshbach resonance 
using various numerical methods for three and four fermions (Stetcu et al. 2007, von 
Stecher & Greene 2007, Blume & Daily 2009, Daily & Blume 2010) and for larger 
systems as well (von Stecher et al. 2008, Zinner 2009). This can be viewed as an 
attempt to address the physics of the BCS-BEC crossover from the few-body side. In 
particular, methods rooted in exact diagonalization that have been used abundantly 
in nuclear physics can be directly applied to the study of fermionic systems as the 
scattering length is varied (Stetcu et al. 2007, Alhassid et al. 2008). 

The shell model Monte Carlo approach (Koonin et al. 1997) is a traditional nuclear 
physics method that has been recently applied to study trapped few-body systems with 
fermions as the interaction is varied (Zinner 2009). In Table [2] we show results for up 
to 21 fermions in a harmonic trap for a positive and a negative scattering length that 
is close to unitarity. Notice how the results for two, three, and four particles lie on 
each side of the energies at unitarity quoted above. The three-dimensional harmonic 
oscillator trap comes with its well-known shell structure. This is illustrated in figure [5] 
where the upper panel shows the energy divided by the scaling expected in a Thomas- 
Fermi approximation. One clearly see an odd-even effect which is most pronounced on 
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Figure 5. Scaled energies Ei n t/N A ^ , where Ei n t = E — E g —o (upper panel), 
and A(N) = E(N + 1) - 2E(N) + E(N - 1) (lower panel) as a function of particle 
number N for a/b = 11 (solid) and a/b = —1.0 (dashed). From (Zinner et al. 2009) 

the BEC side of the resonance. A useful observable for such effects is the pairing gap 
discussed at the end of Section |3.2| In the lower panel of figure [5] the gap is depicted 
as function of particle number. The shells are seen to matter most on the weaker 
BCS side of the resonance. This is consistent with the finding using other methods 
such as fixed-node diffusion (Blume et al. 2007) and Green's function Monte Carlo 
(Chang & Bertsch 2007), and with density functional theory using the local density 
approximation (Bulgac 2007). Similar gaps are also found in two-dimensional systems 
(Rontani et al. 2009). 

We now consider some analogous systems from few-body nuclear physics. 
Nucleons are fcrmionic spin one-half particles and one could imagine having a two- 
neutron or two-proton state, the former would be more desirable due to its lack of 
Coulomb repulsion. However, it is well-known that only the neutron-proton system 
is bound (naturally occuring as the atomic deuteron, 2 D). Considering then three- 
nucleon structures the possiblities are three neutrons, 3 H or 3 He. We are interested 
in a comparison to the few-body cold atom systems discussed above, and will focus 
on neutrons for the moment since the scattering length for neutrons is large. 

Since the neutron-neutron scattering length is so large it indicates that the system 
is close to the resonance (\a\ ~ > oo) where a bound state appears. It is then reasonable 
to explore whether the additional neutrons in three and four neutron clusters could 
provide an overall bound system. Early experimental studies provide no evidence for 
either system (Fairman & Meyerhof 1973). A claim of a tetra-neutron signal appeared 
in 2002 (Marques et al. 2002), which could, however, not be reproduced (Aleksandrov 
et al. 2005). From a theoretical point of view, it was realized much earlier that tri- 
and tetra-neutron cluster are unbound when using nuclear forces that give reasonable 
predictions for nearby bound nuclei such as 3 H, 3 He and 4 He (Bevelacqua 1980). 
Recent calculations using state-of-the-art modern nucleon-nucleon potentials support 
this conclusion for both the tri- (Hemmdan et al. 2002, Lazauskas & Carbonell 2005) 
and tetra-neutron system (Pieper 2003). 

The situation is somewhat different if one starts to consider neutron-rich few- 
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body nuclei that have one or a couple of protons. A nice example is 5 H. It could be 
the structure of a deuteron (neutron-proton) and a tri-neutron but in fact has the 
structure of triton, 3 H, plus two neutrons (de Diego et al. 2007). The latter implies 
that it could have an interesting three-body spectrum related to Efimov states to be 
discussed in Section [6] Adding one or two more neutrons would give hope of realizing 
a tri- or tetra-neutron plus 3 H system. However, these attempts have also been futile. 
In the case of two protons, i.e. in the Helium isotopic chain, one can play the same 
game and this has been continued all the way up to 10 He. A review concerning the 
limits of stability for different isotopic chains can be found in (Thoennessen 2004). 

The nuclear studies of such exotic few-body states with large imbalance of 
neutrons over protons are difficult due to the complicated structure of the nuclear 
interaction. As the interactions are well under control in cold atoms in general and 
also in experiments with only a few fermionic atoms, one can hope to gain insights 
into the general issue of small interacting Fermi systems that can be transfered into 
the nuclear physics domain. 

5. Bosonic structures in few-body systems 

The BEC structure is defined by mean-field properties. The mean-field approximation 
for self-bound finite systems violates translational invariance and includes a spurious 
center-of-mass motion. This approximation easily obscures all finite number effects, 
as e.g. shell effects. Conversely, a proper ./V-body wave function $ only depends 
on N — 1 independent relative coordinates in conflict with the mean-field product 
structure. Suitable coordinates in the one-body density matrix have to be specified to 
allow computation of the eigenvalues which in turn decide whether a given structure 
can be classified as a BEC. We believe that only one combination of coordinates and 
wave functions is reasonable, i.e. use of initial single-particle coordinates and a total 
wave function where the relative wave function \& is multiplied by a center-of-mass 
function adjusted to maximize the largest eigenvalue of the resulting one-body density 
matrix. 

Before application of BEC definitions it is desirable to know whether a given 
system is strongly or weakly correlated as in solids or liquids, crystal or mean- 
field structures, localized or delocalized structures. Second, constraints from basic 
symmetry requirements related to conserved quantum numbers are crucial for finite 
systems. After discussions of these general questions we turn to applications of a- 
clustering and condensation in nuclei. 

5.1. Mean-field versus Localization 

A system of N particles prefer crystallization at low temperature (T) if the two-body 
interaction has a sufficiently deep minimum to localize the relative wave function 
around the corresponding distance. Alternatively, the zero point motion in a weaker 
attraction could be larger than the range of the potential and the particles would 
move across the potential wells created by neighboring particles. Then each particle 
would rather be subject to the average potential from all particles and the mean-field 
approximation would be appropriate. 

To classify the potentials Mottelson, following earlier work by de Boer on the 
solid state of noble gases (Ashcroft & Mermin 1976), introduced in (Mottelson 1999) 
the quantality parameter Ajv/ot as the ratio between the kinetic and potential energy 
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Figure 6. Illustrative plot of the inter-particle potential as a function of relative 
distance. 



of particles located in a relative two-body potential (see Fig. [6]) . Using the notation 
of the figure, this parameter can be written 



where m is the mass of the particles, with Vq the depth and c m i n the position of the 
assumed potential minimum. The kinetic energy K is estimated as K = f? jm,(? min 
with the momentum p ~ ?i/c m ,- n and the reduced mass to/2. The potential energy is 
estimated as the value Vo at the equilibrium position. The quantality parameter does 
not distinguish between bosons and fermions, implying that the criterion is a truly 
profound way to classify physical systems. 

Mottelson noticed (Mottelson 1999) that Amo* ~ 0.1 marks a transition between 
solid and liquid (mean-field) behavior. At T = both atomic 4 He and 3 He are liquids 
with Ajv/ot > 0.1 whereas Ne and molecular H 2 are solids with Am * < 0.1. For 
nucleons A^Mct ~ 0.4 and therefore nuclei are expected to be described rather well 
by mean-field models. The physical interpretation is straightforward. When A^/ ot is 
small, the kinetic energy is small compared to the potential which keeps the particles 
in their equilibrium positions and produces a solid. When A^/ot is large, the kinetic 
energy (zero point energy in the well) is larger than the depth and the particles 
move to distances greater than c m i n . We note that the usual characteristic of a solid 
in condensed-matter physics which is the breaking of translational symmetry in the 
crystal lattice is not obviously applicable to the finite systems we consider here. The 
total wavefunction must of course be rotationally invariant, but by looking at the 
wavefunction as a function of the spatial coordinate for fixed center-of-mass position 
the localized 'solid' structure will emerge as concentrations of the probability density. 
Fig. Q illustrates this for the N = 3 case. 

When we say that the system exhibits mean-field behavior the statement is that 
there are single-particle wave functions extending over the whole volume occupied by 
the particles. However, since we are addressing interacting systems these mean-field 
wave functions include effects from all the other particles. The corresponding particles 
are "effective" in contrast to physical and the wave functions are often called quasi- 
particle orbits. The statement that Ajvfot > 0.1 can then be reformulated by saying 
the low-energy quasi-particle excitations have large (compared to the systems spatial 
extension) mean- free path (Mottelson 1999), which is the same as total delocalization. 



K 



(32) 
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The quantality condition can be related to other criteria from quantum and 
statistical mechanics. We first distinguish between systems at a given temperature 
determined by an external reservoir and finite (nuclear) systems where an unambigous 
way to assign a temperature does not exist. For a specified T the thermal de Broglie 
wavelength for ideal gases is given by (Fetter & Walecka 1971, Huang 1987) 

X dB (T) = J^— , (33) 

which may be obtained from A = 2ith/p and p 2 — 27rm/csT, where p is the momentum, 
m is the particle mass and ks is Boltzmann's constant. This wavelength can be used 
to characterize when a quantum gas of cither fermions or bosons will display strong 
quantum behavior as the temperature is lowered. For a Bose gas the condition is 
1 *C [Ads(7 1 )] 3 n, where n is the particle density. We have d ~ n -1 / 3 (with d the 
inter-particle distance) , and this criterion can thus be written 1 <C XdB (T) / d. For a 
Fermi gas one has similarly XdB/d oc \f up /ksT , where ep is the Fermi energy. We 
thus see that for ksT <C ep, the de Broglie wavelength will be much greater than the 
average inter-particle spacing. 

We used that the kinetic energy of each particle on average is irksT up to factors 
of order unity. We assume that the total energy E of the two-body system is close to 
zero, implying that an estimate of the kinetic energy is \Vo\, see Fig. Q. Then the 
criterion for quantum degenerate behavior is 

1 « A des = ^fl = - 27T \ = y^W, (34) 

where we used d ~ c m i n and V — Vq (see Fig. dfjj)). We therefore see how one 
can relate the quantality condition to the relation for the breakdown of the classical 
gas regime in statistical mechanics. Again there is no distinction between bosons 
and fermions, explaining the universality of Mottelson's criterion for such seemingly 
different systems. 

Alternatively, in a finite isolated system where temperature is not assigned the 
de Broglie wavelength can instead be defined for the motion of the constituents as 
X c ib = 2nh/ ' \]2m(E — V). If we now use E ~ and V — —Vq, we get the same 
relation as in Eq. (34) except for an unimportant factor of y/n. Thus the quantality 
parameter Aj\,/ ot is again related to the de Broglie wavelength of motion, XdB, an d the 
typical inter-particle distance. 

The physics interpretation is that the wavelength, both with and without 
temperature, compared to the inter-particle distance determines the character of the 
structure, i.e. small wavelengths prefer solids and large wavelengths prefer mean- 
field structures. More quantitatively if the wavelength reaches the order of the 
inter-particle distance, coherent quantum structure is preferred. The critical value is 
Adeg — 1 or equivalently Ajv/ot — 1/(2"") = 0.16 in nice agreement with the observation 
that Amo* — 0.1 is the critical value (Mottelson 1999). The competition of mean- 
field and cluster states in nuclei has recently been addressed using state-of-the-art 
nuclear energy density functionals (Ebran et al. 2012a, Ebran et al. 2012&) using both 
relativistic and non-relativistic functionals. The results obtained are consistent with 
the simpler models discussed above. 
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Figure 7. Illustrat ive plot of the delocalized wave function in Eq. 
localized one in Eq. (|36f for N = 3. 
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5.2. Symmetries, Localization versus Derealization 

Precise condensate criteria are derived from the one-body density matrix, pi{f,r l ), 
denned in Eq. ([!]). To illustrate effects of symmetry and localization we use gaussian 
wave functions. This allows analytic derivations yet containing the generic features of 
interest. A perfect mean- field product TV-body wave function is given by 



N 

= (6^)- 3JV/2 exp(-^r l 2 /(26 2 )) , 
1 



(35) 



where ri is the ith coordinate. This wave function, illustrated in Fig. ([7]), as well 
as all other mean-field wave functions violate translation invariance, or equivalently 
momentum conservation, which is restored by integrating ^({r, — R'}) exp(iP ■ R) 
over all R' (Siemens & Jensen 1987). The result, #i n t, of lowest energy has P = 
which for Eq. ([35]) results in Eq. Q which also is invariant under rotations around 
the center-of-mass R. 



In contrast to the mean-field solution in Eq. (35), the particles can be correlated 
and the wave function ^i oc in a body-fixed coordinate system localized in distributions 
around preferred points R^, i.e. 



N 

£ 

i=l 



(r i -R p(t) ) 2 /(2B 2 )), 



(36) 



; ({ r J) « £exp(- 
M 

where the normalization is omitted and full symmetry is achieved by the summation 
over all permutations p of the set of numbers {1,2,...., N}. This is illustrated in Fig. ^ 
for N = 3. The translational invariance is restored precisely as for Eq. (351, i.e. the 
invariant wave function is obtained from Eq. (l36|) by the substitution r j — > qi = r j — R 
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with a corresponding change of normalization constant. The rotational invariance is 



broken for ^i oc in Eq. (36) but recovered for states of zero angular momentum by 



equally weighted linear combinations of all spatial rotations of ^i oc . 

The condensate fraction depends strongly on the degree of localization as we can 



see explicitly by computing the one-body density matrix for Eq. ( 36 1 . We assume very 
narrow non-overlapping gaussians and obtain 

Pl (r,r>) = (BV^r 3/2 (37) 

N 

x Y, exp(-((r - R k f + (r' - R k f)/{2B 2 )) , 

k=l 

which has N equally large eigenvalues while all others are zero. This is a condensate 
fraction of 1/N corresponding to one single-particle state for each of the N particles. 
However, after restoration of rotational symmetry only eigenvalues zero remain. If the 
widths, B, of the gaussians increase and they begin to overlap with each other one 
eigenvalue separates out and becomes finite. Increasing the width leads to increasing 



overlap with a product wave function like Eq. ( 35 ) 



We can quantify by computing the overlap between the factorized and localized 



wave functions in Eqs. (35) and (36), i.e 



(*|*; oc ) = 



2bB 
b 2 + B< 



3N/2 r2 

expf- ^ k=1 k ) 



(38) 



which only is close to unity when b ~ B and either Ri/B <C 1 or R/b <C 1. Eq. (38) 
is also obtained if one averages ^i oc over all possible angles to obtain a rotationally 
invariant wave function. Thus a substantial condensate fraction requires that the 
overlap of ^i oc with Eq. (35) is large. 



The discussion above is very simple to relate to the Mottelson quantality. The 
widths of the gaussians increase with the zero-point energy in the two-body potentials, 
which corresponds to increasing kinetic energy in Eq. ( 32 ) . The states then overlap 



and one eigenvalue of the density matrix will grow in magnitude. On the other hand, 
the depth of the potential tends to localize the particles with decreasing widths for 
increasingly attractive potentials. The density matrix will thus have N eigenvalues of 
order one. 



5.3. Trapped Few-Boson Systems 

As discussed in Section |2.1[ the characterization of cold atomic Bose gases from the 
one-body density matrix has been explored in great detail for large samples of atoms. 
To compare nuclei with relatively few particles to atomic gases we now discuss the case 
of a small number of bosons in a trap. Monte Carlo studies exist for the homogeneous 
(Giorgini et al. 1999) and the trapped case (duBois & Glyde 2001) that consider the 
total energy per particle. However, these results are for large systems of more than 
a hundred bosons. Diffusion Monte Carlo was performed for systems of N = 2 — 50 
particles (Blume & Greene 2001), and also variational Monte Carlo has been used 
(Hanna & Bluhmc 2006). The energy per particle obtained in the smaller systems are 
in good agreement with that obtained from larger systems. However, they also show 
considerable depletion of the condensate fraction at large scattering lengths. Common 
to these studies is the use of repulsive two-body potentials. This has the drawback that 
the range of the potential increases essentially linearly with the scattering length, i.e. 
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at large scattering length the range of the potential can become as large as the inter- 
particle spacing and the details of the potential can no longer be neglected. In this 
non- universal limit, where the range of the potential is comparable to the scattering 
length, the energy cannot be expressed in terms of only the scattering length, and an 
effective one-channel model for the effective interaction is no longer sufficient. 

A solution to this problem is to use an attractive finite-range interaction as 
a model potential. This would also seem more realistic in terms of providing an 
effective one-channel model for the Feshbach resonances of experiments. It is well- 
known that potentials like the attractive gaussian or square well can attain any value 
for the scattering length as one increases the strength around the threshold for the 
appearance of a two-body bound state. This brings in the complication that the many- 
body problem will have many self-bound negative-energy states when bound states 
are allowed for a pair of bosons. However, if one uses direct numerical diagonalization, 
then all states are determined and it turns out that states with the properties expected 
of condensates appear as excited states in a quasi-continuum that are orthogonal to 
all the bound negative-energy states (Th0gersen et al. 2007). An example of this is 
shown in Fig. [S] 
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Figure 8. A typical spectrum (in the vicinity of zero energy) of a system 
of N = 20 bosons in an isotropic harmonic oscillator trap, interacting via an 
attractive gaussian two-body potential with one bound state and a positive 
scattering length. The inset shows the beginning of the so-called quasi-continuum 
spectrum. From (Th0gersen et al. 2007). 

The condensate nature of the lowest state in the quasi-continuum can be verified 
by calculating the central density of the cloud. For a completely delocalized mean-field 
condensate state one would expect the bosons to occupy the entire trap. Indeed, as 
shown in Fig. [9] one finds that the first state of positive energy fulfills this conditions, 
whereas the negative-energy self-bound states have a much larger density and thus 
smaller volume per particle. This is also reflected in the condensate fraction which 
is small for the negative-energy states, is unity for the lowest state in the quasi- 
continuum, and then decreses slowly for higher- lying states (Th0gersen et al. 2007). 

Another very important question is the behavior of the condensate fraction as 
a function of particle number. In Fig. [10] the condensate fraction is shown as the 
(positive) scattering length is increased for 10, 20, and 30 bosons in a trap. We clearly 
see the depletion of the condensate as the interaction is increased, and that this effect 
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Figure 9. The inverse central density n 1 (in oscillator units n OBC = 
n~ 3 / 2 Nb~ 3 , where b is the oscillator length) for a system of N bosons in an 
isotropic harmonic oscillator trap, interacting via an attractive gaussian two-body 
potential, with scattering length a = 119au as function of the state number. The 
lowest state with positive energy is numbered zero. From (Th0gersen et al. 2007). 



is more severe as the particle number decreases. This demonstrates the general trend 
that small condensates are very fragile when interactions, and in turn correlations, are 
strong. Similar results can be obtained using harmonic approximation methods to the 
TV-body problem in one, two and three dimensions at the Hamiltonian level (Zaluska- 
Kotur et al. 2000, Yan 2003, Gajda 2006, Armstrong et al. 2011, Armstrong et al. 
2012a, Armstrong et al. 20126, Armstrong et al. 2012c, Volosniev et al. 2013) and via 
path integrals (Brosens et al. 1997a, Brosens et al. 19976, Brosens et al. 1997c, Brosens 
et al. 1998a, Brosens et al. 19986, Tempere et al. 1998, Tempere et al. 2000, Lemmens 
et al. 1999). These facts should serve as a warning that systems with very low particle 
numbers are hard to meaningfully describe in the same terms as condensates. Below 
we will look at an example from nuclear physics and the use of a-particles as a degree 
of freedom in the nucleus. 

5.4- Nuclear a -particle condensates 

The idea of a-particles as essential entities in the structure of nuclei arises from the 
small a-radius, the relatively large binding and the spin saturation of both neutrons 
and protons. Attempts were made in the early days of nuclear physics to construct 
nuclear structure from a-particles and valence nucleons (Wigner 1937, Wheeler 1937, 
Wefelmeier 1937). When more information became available these cluster structures 
turned out to be less than robust as indicated by the overlapping densities of the 
different a-particles. Thus there are no compelling reasons for clusterization of 
nucleons into a-particles in the bulk at nuclear densities. A review on clustering 
in general and a-clusters in particular can be found in (Freer 2007). 

On the other hand, a-cluster models are able to explain many properties of 
specific light nuclei (Brink 1966). The most prominent examples are the ground 
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Figure 10. The condensate fraction of a BEC-state of a system of TV identical 
bosons in an isotropoci harmonic oscillator trap with attractive gaussian two-body 
potential as function of (TV — l)(a/fe) 1//2 where b is the oscillator length. From 
(Th0gersen et al. 2007). 

state of 8 Be and the so-called Hoyle-state in 12 C, both a-unstable + -states below 
all other particle emission thresholds (Hoyle 1954, Cook 1957). These structures were 
accurately described many years ago (Horiuchi 1970, Harvey & Jensen 1972, Uegaki 
et al. 1977) and later on confirmed in numerous theoretical works, see e.g. (Kanada- 
En'yo 1998, Descouvemont 2002, Neff & Feldmeier 2004). Based on calculations 
that are approximations to earlier results these states were recently interpreted as 
a condensate of respectively two and three a-particles (Tohsaki et al. 2001, Tohsaki 
et al. 2004, Funaki et al. 2002, Funaki et al. 2008, Funaki et al. 2009, Funaki et al. 2010). 

Since the a-particles have spatial extension and internal structure, the general 
condensate condition must be supplemented by the further requirements (i) that inert 
a-particle constituents must be present, and (ii) that the a-cluster wave function must 
be orthogonal to all other many-body excited states or equivalently that the internal 
a-particle degrees of freedom decouple from relative degrees of freedom. 

The quantality condition for an effective a — a potential can be estimated through 
a common parametrization in terms of an attractive and a repulsive gaussian (Ali & 
Bodmer 1966). The minimum value is Vq ~ 5 — 8 MeV (including the Coulomb energy 
of « 2 MeV) for c min w 2.5 — 3.0 fm and Am ot 0.1 — 0.2. Thus the mean-field model 
is first choice for point-like a-particles but they have a strong tendency to localize. 
Thus a substantial condensate fraction can be expected for such idealized point-like 
particles. 

However, a-particles have internal structure, finite extension, and interact via 
both attractive short-range and repulsive long-range forces. They have essentially 
to touch each other to form a resonance or bound state because the Coulomb 
repulsion otherwise would blow the pieces apart. At such relatively small distances the 
intrinsic nucleonic many-body degrees of freedom easily get excited. The restricted 
Hilbert space describing condensates, i.e. mean-field for a-particles and all nucleonic 
structures frozen, gets increasingly inadequate with increasing excitation energy, E*, 
simply because the many-body density of states increases exponentially. The a-cluster 
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structure is most likely to be prominent at the a-disintegration thresholds which 
increase with nuclear excitation energy as E* « 7(N a — 2) MeV, where N a is the 
number of a-particles in the cluster. Thus only small N a is possible. Quantitative 
estimates can be found in (Zinner & Jensen 2008a). 

To be specific we illustrate with details from the most well-known candidate, 
i.e. the Hoyle state of three a-particles in 12 C which has about 90% overlap with 
a-particles in relative s-waves around the center-of-mass (Matsumura & Suzuki 2004, 
Suzuki & Takahashi 2002, Funaki et al. 2003, Chernykh et al. 2007). This corresponds 
to a one-body density eigenvalue of about 0.7, computed by assuming the optimum 
center-of-mass wave function (Matsumura & Suzuki 2004, Suzuki & Takahashi 2002). 
At the same time a-cluster models show a-particle density distributions localized 
around specific points in space (Chernykh et al. 2007). Reconciling these results, 
where apparently both localization and large condensate fraction are present in the 
same wave function, is only possible with large widths of the localized wave in Eq. ( 36 ) . 
This effectively recovers (most of) the independent particle wave function in Eq. (35). 
Since the width is large the particles are sufficiently separated to remove the need for 
nucleon antisymmetrization and the a-particles are present. 

Thus, condition (i) is rather easily fulfilled provided the wave functions are 
constrained to describe essentially non-overlapping a-clusters. However, this structure 
is rarely an accurate approximation as conditions (ii) only is expected to hold for 
N a < 4 for the condensate candidates. For low excitation energy, the a-particles would 
overlap and the nucleon degrees of freedom wash out the a-clusterization. At higher 
excitation energy, the condensate wave function constructed in the limited Hilbcrt 
space, would be smeared out over numerous other close-lying states. Collective states 
are also very unlikely to display features of a condensate as the Coulomb interaction 
will mix these with other states. The realization of condensates of a particles in nuclei 
therefore face severe and likely insurmountable obstacles. 

In relation to our discussion of trapped few-body Bose systems in Section [5.3[ we 
note that for three or four bosons the depletion for strong interaction is severe. The 
a-particle studies use an effective a-a interaction that is typically a parametrisation 
as in (Ali & Bodmer 1966). Small changes in this interaction can thus cause large 
deviations in the condensate fraction. Describing these a-cluster states in terms of 
properties applied to condensates can therefore quickly be very model-dependent, 
leaving the interpretation less useful. 

The linear chain structures of a-particles (Zinner & Jensen 2008a) at the break- 
up threshold (which can be organized in so-called Ikeda diagrams (Ikeda et al. 1968)) 
are conceptually similar to the one-dimensional atomic condensates called Tonks- 
Giradeau structures (Paredes et al. 2004, Kinoshita et al. 2004). Such nuclear states 
have been searched for and for many years the Hoyle state in 12 C was the favorite 
candidate. However, the linear chain suffers from precisely the same difficulties as the 
a-condensate discussed above (Zinner & Jensen 2008a). 

Discussions about possible condensate structure are not constructive without 
an experimentally observable distinction from other structures. First, all one-body 
properties are excluded as signals. Condensates are diagnosed through coherence 
properties of the wave function, not by density distributions. This immediately 
implies that one-body observables like elastic scattering cross sections used in 
(Chernykh et al. 2007) cannot be used to distinguish between condensates and other 
structures. Second, the experiments almost inevitably at some point detect large- 
distance properties which in the context of coherence is related to the long-range order 
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parameter. For a-particles no product structure at large distance is possible due to the 
Coulomb repulsion. To see the short-distance bulk structure responsible for the large 
one-body density eigenvalue require detection of properties of structures differing by 
one or more particles (Gajda 2006). The results should remain unchanged except for 
the normalization proportional to particle number. Based on these observations, it is 
unclear what new physical consequences arise from approximating the fully correlated 
cluster states (as in (Uegaki et al. 1977)) with less accurate a-particle condensate 
states as no observable are available that can distinguish between the two. In any 
case one must also remember that all cluster states are approximations to the real 
nuclear states. 



6. Efimov Physics 



About forty years ago the Efimov effect was suggested as an anomaly appearing in a 
three-body system when at least two of the two-body subsystems simultaneously have 
short-ranged bound states at zero energy (Efimov 1970, Efimov 1971). The three-body 
system then has infinitely many bound states with exponentially decreasing binding 
energies and correspondingly increasing mean square radii. In figure [TT] we present an 
example of the three-body spectrum considered by Efimov. The properties of these 
states are independent of the potentials supporting the two-body states at zero binding 
energy. The only requirements are that three dimensions are available, the potentials 
are of short-range and sufficiently attractive to provide borderline binding. We shall in 
this section give the simple arguments with corresponding derivations, and generalize 
the concept of model independence to a broader class of phenomena called Efimov 
physics. We shall discuss possible occurrence in nuclear, molecular or atomic systems, 
and the recent experimental verification in cold atoms using the technique of Feshbach 
resonances. 



6.1. Basic derivation and occurrence conditions 

Convenient coordinates for three particles are the hyperradius p and the five 
hyperangles £1. It suffice here to define p by 

mp 2 = m,m,(rj — r, ) 2 (39) 

J i<i 

in terms of an arbitrary normalization mass m, the three masses m^, k = 1,2,3, 
and the three particle coordinates r^. The pertinent properties of the Schrodinger or 
Faddeev equations in the present context can be extracted from the lowest effective 
(hyper)radial equation (Nielsen et al. 2001) 



d^_e ± i/A , _ , 

dp 2 p 



,2 



K 



f(p) = , (40) 



where k 2 = —2mE/h 2 , and £ is a dimensionless function of p determined from the 
interactions and the angular equations. Outside the ranges of the interactions the 
function £ turns out to be independent of p, strongly mass (m^) dependent (Jensen 
& Fedorov 2003), and £ is a real number provided two or all three subsystems have 
bound states at zero energy. 

The solution for the radial wavefunction is then / oc ^fpKi^np) where is the 
modified Bessel function of the second kind. For small k we get / cx v /psin(^ ln(p/po) 
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Figure 11. Energy of Efimov States as function of the scattering length, k* 
is the background scale fixing the Efimov spectrum as explained in the text. On 
the a < side there are no bound dimers and the trimer states (T) terminate in 
the three-particle continuum (A + A + A), whereas on the a > side the states 
terminate at the atom-dimer continuum (A + D). Both axis have been scaled by 
a power of 1/8 in order to make the similiarities in the spectrum clear. 



where po is a constant defining the arbitrary length scale chosen from eq.(40) when 
E — > 0. Then po is the first node of the wavefunction. This app roximation for / could 
be found directly by testing the solution / oc p x l 2±% ^ in eq.(40) for small k. For large 
p the wavefunction must fall off exponentially for a bound state, i.e. / oc exp(— up). 
The bound state energies are found by the condition Ki^(Kpo) = 0, that is from the 
nodes of the Bessel function K^(z n ) = with n n = z n /po. For small k this implies 
z n = z exp(-mr/£) or 

E n = g exp(-2mr/0 , (p 2 ) n = ^ ^ ■ (41) 

These scaling relations clearly only depend on the value of £ and an energy or length 
scale Eq or po which are determined by the finite-range potentials and the boundary 
conditions at small distance. If the potentials are strongly attractive a number of 
bound three-body states may appear at energies lower than the scaled sequence in 
eq.(|4l]). 

This scaling limit is reached when the scattering lengths diverges to infinity 
(bound state at zero energy) in comparison with the ranges of the potential. The 
sequence is terminated when the spatial extension of the state becomes comparable 
to an average of the scattering lengths. Thus we find the Efimov sequence for large 
scattering length, as well as for potential ranges approaching zero which is the Thomas 
collapse (Thomas 1935). Both effects are therefore described by these expressions. 

The crucial property is the inverse square behavior of the effective potential in 
eq.(40) in (at least) a region of space. This form arises at distances between the 
effective range R e of the two-body potentials and their scattering length a. All Efimov 
states are located in this region of space. By counting the number of nodes of / in this 
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interval we find the number of bound states to be Ne ~ C/ 71 " m ( a e/ Re)- If th e strength 
of the 1/p 2 potential is —1/4 only an infinitesimal negative or positive addition would 
produce either infinitely many bound states or none at all. 

The decisive parameter is then the strength £ which determines the character of 
the potential, that is infinitely many for real £ or no bound state for imaginary £. 
The eigenvalue of the hyper-angular equation determines £. For identical bosons the 
equation is (Jensen & Fedorov 2003) 

8sinh(£7r/6) = £\/3 cosh(^7r/2) . (42) 

with the solution £ = 1.00624. For non-identical bosons where all scattering lengths 
still are large the equation is (Jensen et al. 1997) 



£ cosh(^/2) \ 3 /£ cosh(£7r/2) \ (ff + /£ + /|) 



2F J V 2F J F 2 
where F = (/1/2/3) 1 / 3 and 



= 2 , (43) 



sinh(£(7r/2 - tp k )) 

Jk = i—n. x (44) 

sm(2ip k ) 

ifik = arctan ^Jmk(m.i + nij +m( : )/(m 1 mj)j . (45) 



The limit of identical bosons is obtained from eq.(43) for (p k — tt/3. When only the 
two scattering lengths between particle pairs (j, k) and (i, k) are large the £ equation 
becomes (Jensen & Fedorov 2003) 

£ cosh(£7r/2) sin(2<^) = 2 sinh(£(7r/2 - tp k ) , (46) 

where ipf. is now the only mass dependent parameter. The solution for s = £ as 



function of the two mass ratios, mijm\ and m^/rni, is shown in figure 12 

The variation of £ with the mass ratios cover all values from zero to infinity. 
Small values imply large distances between the Efimov states and for realistic finite 
scattering lengths there is only room for very few with smaller binding. For identical 
bosons the ratio of radii between two neigboring states is 22.7, that is to allow three 
states the scattering length must be about 10.000 times larger than the range of the 
potential. 

On the other hand large values of £ is also possible and many Efimov states can fit 
into the allowed interval. This is achieved when two masses are much larger than the 
third. The opposite situation with two light and one heavy mass is the least favorable 
to find Efimov states. Since at least two scattering lengths must be large it is much 
simpler to use two identical particles where any required tuning automatically is done 
simultaneously for both particles. However, this may also reduce the contributions to 



only these two subsystems where eq. ( 46 1 applies instead of eq. ( 42 ) . If we maintain 
equal masses we find the somewhat less favorable value of £ = 0.499 instead of 
£ = 1.00624. However, if the mass of particle k is small compared to the masses 
of the other particles the value of £ can be very large. The mass variation is much 
more efficient than increasing from two to three contributing subsystems. 

The quantum statistics due to identical particles does not change anything for 
bosons. In contrast, none of the derivations are valid for three identical fermions 
where the antisymmetry requires either different s-states or at least one p-state. For 
two identical fermions precisely the same formulations apply since each spatial state 
can hold at least two fermions by use of different spin projections. The antisymmetry 
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Figure 12. Surface plot of the scaling factor, so = £j as function of the two 
mass ratios, mil 'mi and ms/mi. The mass ratios are on a logarithmic scale. 

is taken care of by the spin degree of freedom which effectively is absent for spin- 
independent interactions. 

Unfortunately the jargon in present days cold atom physics can easily cause some 
confusion as we have discuss above. The single-particle states are hyperfine states 
of given angular momentum projection and a set of other quantum numbers related 
to the trap or momentum vector for a homogeneous system. The spin projection 
is then replaced by the hyperfine quantum numbers. Thus, since antisymmetry still 
is required, s-wave contact interactions are forbidden between particles in the same 
hyperfine state. Since two subsystems in s-waves are necessary with large scattering 
lengths we must have at least two different hyperfine states. Thus fermions in one 
of these states and the third particle (boson or fcrmion) in another state leave one 
p-wave interaction and two possible s-wave contact interactions. The total angular 
momentum is then at least 1 unit of h. It is of course also possible to have identical 
fermions in three different hyperfine states. This allows all pairs to interact by s-wave 
contact interactions (Williams et al. 2009, Wenz et al. 2009, Nakajima et al. 2011). 

The question then arises if Efimov states exist for finite three-body angular 
momentum L > 0. The answer is yes (Nielsen et al. 2001) but with conditions of 
sufficiently small mass, Wfc, for the third s-wave interacting particle. The critical 
mass depends on L, which is the same as the orbital angular momentum of the two 
non-s wave interacting particles. For L = 1 the result was given as 

mfc < 0.5(— rrii — rrij + <J (mj + raj) 2 + 0.60mirrij . (47) 

Again only mass ratios enter the condition. This result is valid independent of the 
boson or fermion character of the particles. In particular, with all masses equal the 
requirement cannot be fulfilled since 1 > 0.5(— 2 + \^4~6. Thus equal masses cannot 
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produce L = 1 Efimov states. Other studies have addressed the case of three identical 
fermions which require higher-partial wave interactions to accomodate the necessary 
antisymmetry (Suno et al. 2003, Macek & Sternberg 2006, Jona-Lasinio et al. 2008). 
We note that very recent studies indicate that p-wave interactions produce unphysical 
negative probabalities at the two-body level and thus also rule out an Efimov effect 
with p-wave interactions (Nishida 2012, Braaten et al. 2012). 

As a curiosum we mention that the formulations of occurrence conditions can be 
investigated in any dimension d. The mathematical equations only leads to potentials 
of form eq.(40) and real values of £ when 2.3 < d < 3.8. Thus, d — 3 is the only 
integer dimension to exhibit the Efimov effect, and for that matter also the Thomas 
effect. However, recent investigations softened these conditions by allowing different 
participating particles to move in different dimensions, e.g. two particles confined to 
d = 2 and a third in d — 3 (see (Nishida & Tan 2011) for a recent review of such 
possibilities). 

In a strict two dimensional setup we do not have the Efimov anomaly with 
infinitely many bound states around threshold but the spectrum of three-body states in 
the zero-range limit still has some rather peculiar features (Bruch & Tjon 1979, Nielsen 
et al. 1997, Nielsen et al. 1999, Volosniev et al. 2012a). In fact, for a system of three 
identical bosons, one finds exactly two bound three-body states, a ground state with 
energy E 3 = 16.52E 2 and an excited state with energy E s = 1.27E 2 . Here E 2 < 
is the energy of the two-body bound state of the zero-range potentials which always 
exists (Landau & Lifshitz 1981). Bosonic few-body states in two-dimensional setups 
with three (Adhikari et al. 1988, Adhikari et al. 1993, Brodsky et al. 2006, Kartavtsev 
& Malykh 2006) and four particles (Brodsky et al. 2006, Platter et al. 2004), as well 
as larger cluster have been studied (Hammer & Son 2004, Blume 2005, Lee 2006). 
In addition, the scattering (Brodsky et al. 2006) and recombination (Helfrich & 
Hammer 2011) has also been discussed. Recently, it explored how the number of bound 
states changes in system of two identical particles and a third distinct particle and for 
three distinct particles (Bellotti et al. 2011, Bellotti et al. 2012, Bellotti et al. 2013a). 



6.2. Scaling properties and examples of Efimov physics 

The almost complete absence of specific potential dependence of the derived results 
strongly suggests that all observables are connected by universal relations. Such model 
independent results are obviously valuable as predictions virtually without any input, 
apart from validity conditions, and as reference units to evaluate similarities and 
differences. A prominent and relevant example in the present context is early nucleon- 
nucleon scattering experiments where disparate potentials were able to reproduce the 
measured low-energy cross sections. The explanation was given by the effective range 
expansion theory where the phase shift 6 was expanded as function of energy, i.e. 

fccottf(fe) = -- + \k 2 R{k) , (48) 
a 2 

and E = Y?k 2 j (2[T). The scattering length, a, is a constant and R approaches a 
constant, the effective range, r e0 , for small k. Since a is a quantity depending on 
the entire potential similar to a weighted integral, the same low-energy cross section 
proportional to a 2 can be found from tremendously different potentials. Thus, all 
low-energy observables can be expressed by a without reference to specific potentials. 
It should be emphasized that the scattering length is periodic as function of potential 
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strength. The same a-dependence emerges independent of the number of deeper lying 
bound states. 

A related result is the mean square radius of a weakly bound two-body state, 
i.e. < r 2 >— h 2 / (Ap\E\), which only depends on the energy and the reduced mass p. 
This radius-energy relation suggests a similar correlation for three-body systems. The 
mean square radius is now uniquely defined by < p 2 >— N < r 2 > but now far from 
being a uniques function of the three-body energy. First of all the quantities are not 
dimensionless but this is a problem already for the two-body system. 

The first attempt towards universal curves is to use an appropriate length R 
related to the range of the potentials and consider < (r/R) 2 >= ti 2 / (4p\E\R 2 ). Then 
all information about both radii and energies are maintained and the scale of the 
system does not enter. For short-range potentials a square well radius equivalent to 
the range of the true potential would be tempting. Since the spatial extension usually 
is measured by the root mean square value we could choose identical second moments 
of the potentials as criteron which for a gaussian of range b implies that R = b-Jb/2. 

For three-body systems we still have to choose the scale po (see below for some 
recent developments regarding this cut-off). One simple option would be to use 
eq.(39) with p — > p and the two-body square well equivalent radii substituted on 
the right hand side. The curves of < [p/po) 2 ) > as function p\E\pj j /h 2 is now 
closer to being universal. However, two different universal scaling functions appear 
(Jensen et al. 2004, Yamashita et al. 2011). One corresponds to the Efimov states with 
radii and energies related through eq.(41) either by varying £ or for fixed £ following 
the sequence of excited states. This is for large \a\. The other universal function 
corresponds to size and energy of Borromean states approaching the threshold for 
binding which occurs for moderate values of \a\. 

For non-identical particles the individual two-body interactions can be varied 
independently to approach different thresholds for binding. Completely different 
behaviors, although still universal, are then possible for example transitions to two- 
body scaling behavior must appear in the limit towards a bound two-body subsystem. 

Another established universal curve is the Phillips line, i.e. the three-body binding 
(triton) energy as function of the scattering length of the particle-dimer (nucleon- 
deuteron) system. This universality can be shown to originate from small energies 
of all involved subsystems in comparison to the nuclear potentials. The advantage 
of both radius and energy information can be traded for a proper universal curve. 
The ratio of two consecutive three-body binding energies (|-E' n+ i|/|£ , n |) as function of 
relative two to three-body energies, B-2/\E n \, was suggested in (Frederico et al. 1999). 
Measured points on such a curve are then signals of structure similar to the Efimov 
states. 

The Efimov effect is concerned with the anomaly found for three particles when 
two-body subsystems are at the threshold of binding. The results are universal or 
model independent by depending only on the masses. This window of universality is 
only open when the scattering length, \a\, is large compared to the effective range of 
the interaction. It is then interesting to extend to other systems with a similar model 
independent description. The broader class of phenomena is naturally named Efimov 
physics which could be defined as physics where Universality and Scale Invariance 
apply. To be specific Universality then means that one global parameter determines 
all properties of the TV-body system. Scale Invariance means that the same properties 
appear for any length scale whether it is in nuclei, atoms or molecules. 
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These concepts are easily misinterpreted to interchange meaning, i.e. scale 
invariance is taken to mean that the next state in a Efimov sequence of a given 
system has the same property as the previous state, and universality is taken to mean 
that different systems can be described with the same theory. We shall use our first 
definitions where universality means model independence and scale invariance means 
that the theory applies to any length and energy scale. The two ingredients begin to 
be similar when dimensionless quantities are used to express the physics observables. 
The scale dependence disappears and universal features are isolated. 

The window of universality is for three particles between the two-body effective 
range and the scattering length. In terms of the hyperradius we have r e < p < a. 
Viewed from the adiabatic hyper-radial potential the states are in the attractive 
region with positive hyper-radial kinetic energy However, the size of the system 
is large and the two-body subsystems are quickly found far outside the range of their 
attraction. This means that classically forbidden regions are occupied with negative 
two-body kinetic energy. The origin of the model-independent features are then simply 
understood as properties determined outside the region of non-zero potential where the 
kinetic energy term decides the behavior. Due to quantum mechanics, the boundary 
conditions still link to the wavefunction behavior at small disance, but this is then the 
entire dependence describable by the model-independent scattering length parameter. 
Successful descriptions require spatially extended cluster systems which is difficult 
to achieve for self-bound systems with short-range interactions. Only s- and p-wave 
systems allow this. The other direction of more detailed information require short- 
distance properties which essentially excludes zero-range models even when effective 
range improvements are employed. 

A number of systems are already found to exhibit universal features. The four- 
body energy (a-particle) is a unique function (Tjon line) of the three-body energy, 
i.e. the three and four-body energies are always found on the same curve (Tjon 1975). 
This result emerge from zero-range models which must be regularized to avoid Thomas 
collapse. It was shown in (Fcdorov & Jensen 2001) that stabilization of the three-body 
system automatically ensures finite results for the iV-body system. However, this does 
not immediately imply that all TV-body systems are determined from two and three- 
body properties. This is still an open problem (Yamashita et al. 2010). 

Universal behavior of excited four-boson states was recently found in both 
variational (von Stecher et al. 2009) and zero-range models (Hammer & Platter 2007) . 
Each three-body state has two related four-body states at larger binding energy than 
this three-body state (Hanna & Bluhme 2006, von Stecher et al. 2009). They have 
too little energy to decay into this three-body state. In contrast, it was found in 
(Yamashita et al. 2006) that the four-body state depends on interaction details. 

The concept of halo nuclei has been rudimentary extended to more than three 
particles. In general fully model independent nuclear TV-body structures for N > 3 
cannot be expected since all clusters are charged and either confined to small distances 
or repel each other by the Coulomb interaction. The line of arguments begins 
with assumptions of only s-waves and no correlations which prohibits halo existence 
(Riisager et al. 2000, Jensen et al. 2004). If correlations develop the iV-body system 
tends to form substructures and effective reduction of the number of clusters. It is 
therefore interesting that (Yamashita et al. 2010) found that four and five identical 
bosons converge to a radius larger than the interaction range as the threshold of 
binding is approached. These weakly bound states have universal structures. More 
generally, it may be possible to use Efimov-type structures for N = 3 to build general 
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so-called higher-order Brunnian states (Baas et al. 2012). 

Another class of universal Efimov-like structures was found in (S0rensen et al. 
2002, Th0gersen et al. 2008a, von Stecher 2010) for TV-body systems with the 
assumption of only two-body correlations in the variational wavefunction. This does 
not contradict the conclusion in (Amado & Noble 1971, Amado & Noble 1973) that 
the Efimov effect does not exist for N > 3. This result is more specifically stating 
that the threshold for binding the TV-body system is not an accumulation point for 
infinitely many states. 

For dilute systems two-body correlations are expected to be dominating. The 
sequence of excited states found between the interaction range and the scattering 
length obey the Efimov scaling equations in eq.(41|. Whether these states 
(approximately) remain after extension of the Hilbcrt space to allow all correlations 
still needs to be tested. Where this conservation of identity is most likely to be 
preserved, for large scattering length or close to threshold of binding the A-body 
system, is unclear at the moment. Ground and lowest excited states may be outside 
the universal region and even higher excited states are the only hope for finding 
universal behavior. The scaling properties are expected to differ in these regions, 
as e.g. Borromean three-body states scale differently from Efimov states. In any 
case, Coulomb interactions between the particles strongly suppress the possibilities, 
see Section [5] 



6.3. Measurable consequences in physics systems 

In nuclei the clusters are all charged except for systems with neutrons surrounding 
a core. The Coulomb interaction is long-ranged and prevent the Efimov effect as it 
reaches to larger distances than the 1/p 2 potential. In the pure form we are left with 
two neutrons and a core which again is unfavorable since the core necessarily is heavier 
than the neutrons. The well known case is u Li with approximately three scattering 
lengths each of about 20 fm which is about 8 times the core radius. The £ value is 
0.074 and the scale factor on energy and square radius is about 3 ■ 10 18 . One Efimov 
state could then by chance be seen in the interval, but not two, and the scaling can 
thus not be tested. 

Realistic two-body interactions for the two neutron- 9 Li system produce one bound 
state which moves down to about 1 MeV below threshold when the neutron- 9 Li 
interaction is tuned to have a bound state at zero energy. However, the many excited 
Efimov states can not be seen since already the first is located at an enormous distance 
outside the radius of the nucleus (Fedorov et al. 1994). This calculation assumes that 
9 Li has angular momentum zero and therefore precisely the same interaction with 
both neutrons. For a non-zero 9 Li-spin the hyperfine splitting of the strong interaction 
prevents simultaneous large scattering lengths for both neutrons. Thus the Efimov 
conditions are therefore most easily met in nuclei for two neutrons and an even-even 
core-nucleus with zero spin. 

It is still possible to see reminiscenses of the Efimov effect in nuclear decays 
and possibly in nuclear scattering (Garrido et al. 2006). Assume for example 
that 11 Li is excited to a T state above threshold for two- neutron emission. The 
momentum distributions of the decay products, two neutrons and 9 Li, can be measured 
simultaneously. Then characteristic peaks should appear in the probability for 
emission of both low and high energy neutrons. This corresponds to one neutron 
correlated with 9 Li at high energy and the other then emerges with very low 
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energy. The energy distribution for 9 Li has a peak at high energy corresponding to 
correlated emission of the two neutrons and a peak at intermediate energy originating 
from correlated emission of the neutron- 9 Li system. These structures are found 
experimentally (Hagino et al. 2009). 

An example of an unexpected consequence is in dense plasma of a-particles and 
electrons (Jensen et al. 1995). The electron screening reduces the Coulomb repulsion 
between the a-particles, and in particular the long-range Coulomb tail is removed. 
This has first two interesting consequences, i.e. the energy of the 8 Be nucleus decreases 
towards zero and the Coulomb tail decreases between the a-particles. This in turn 
implies that the Efimov conditions are approached, and the corresponding three-body 
system, 12 C, should have decreasing energies and if the screening is sufficiently strong 
more states should appear at the threshold. This is the triple a process proceeding 
via the Hoyle state. The rate could easily be dramatically changed for high plasma 
density and temperature. 

An even more speculative example is the suggestion of recombination of three 
deuterons where two deuterons are bound in a molecule in a dense lattice (Engvild 
& Kowalski 2006). The third deuteron is injected and, catalysis- like, via three-body 
long-range interactions cause fusion of the molecular deuterons and expulsion of the 
injected deuteron. The energy gain is from formation of the a-particle. 

The possibilities are enhanced by use of asymmetric systems with two heavy 
particles. Three-body combinations of 3 He-atoms and two alkali atoms could 
potentially increase the ^-values substantially. As the scaling is exponential a number 
of Efimov states could appear with energy ratios of only 5 — 10. The mass ratio can 
be made much more extreme by using an electron as one of the particles and atoms 
or molecules for the other two (Jensen & Fedorov 2003). Then the energy and radius 
ratio of consecutive states could be very much more favorable by differing only a few 
percent. 

Efimov scaling related to 1 /r 2 potentials appear obviously already for a charged 
particle in a sufficiently strong dipole field. However, coupling to rotational degrees of 
freedom would quickly destroy the simple properties. It was recently suggested that 
quantum dots and artificial atoms could prevent this rotational coupling and thus 
maintain the Efimov sequence (Schumayer et al. 2010). 

Within the last five years a number of papers reported on experimental results 
obtained for cold atoms in traps and with Feshbach resonance techniques (Kraemer 
et al. 2006, Ottenstein et al. 2008, Zaccanti et al. 2009, Gross et al. 2009, Pollack, 
Dries & Hulet 2009, Barontini et al. 2009, Knoop et al. 2009, Huckans et al. 2009). 
The systems always consist of many particles and the Efimov states are only indirectly 
observed as enhanced (or depleted) probabilities for three-body decay of the trapped 
particles. The two-body effective interaction is varied through the magnetic field to 
pass values where the three-body system has zero energy. The decay into a deep dimer 
and a high energy third particle is enhanced since the relative and total energy of all 
three particles in the Efimov state is identical to their energies in the ultra cold trap. 
The coupling between these states is maximized and there is a transition through the 
Efimov state at threshold into deeper-lying dimer states accompanied by the third 
particle with the surplus energy to ensure energy conservation. 

For very attractive two-body interactions that have positive scattering lengths 
(a > 0) and where a high-lying or shallow so-called Feshbach dimer state is present, 
a somewhat similar process also is enhanced. Two states have the same energy, that 
is the well-bound three-body Efimov state and a bound dimer with zero binding 
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of the third atom. From an atom-molecular cold gas the coupling between these 
states is largest and decay through this Efimov state into deeper-lying dimers is 
maximized. Periodic loss rate minima were predicted (Fedichev et al. 1996, Esry 
et al. 1999, Nielsen & Macek 1999, Nielsen et al. 2002) and recently also observed 
for two neighboring structures (Zaccanti et al. 2009). These variations should 
appear periodically with sizes scaled by the Efimov scaling factor of 22.7 for equal 
mass bosons. This three-body recombination process is discussed in details in 
(Esry & D'Incao 2007) and subsequently applied in analysis of the experiments. A 
number of these features are now confirmed (Pollack, Dries & Hulet 2009, Zaccanti 
et al. 2009, Gross et al. 2009), including the scale factor. We note that it is now 
also possible to study Efimov states in three-component 6 Li gases via radio-frequency 
association (Lompc et al. 2010, Nakajima et al. 2011). 

The experimental study using 39 K presented in (Zaccanti et al. 2009) reported 
some deviations from universality which could potentially come from finite-range 
corrections beyond the universal theories that only take the scattering length into 
account. Some recent theoretical works have addressed the corrections coming 
from the finite range of inter-atomic potentials in the coordinate-space formalism 
(Th0gersen et al. 20086, Th0gersen, Fedorov, Jensen, Esry & Wang 2009, Wang, 
D'Incao & Esry 2011, S0rensen et al. 2011) and using momentum-space effective field 
theory (Massignan k Stoof 2008, Platter et al. 2009, Ji et al. 2010, Naidon et al. 2012). 
In particular, the study of (Massignan & Stoof 2008) provide results that are close 
to the experimental data. On the nuclear physics side, finite-range corrections for 
three-nucleon systems have been considered by Efimov himself about two decades ago 
(Efimov 1991). 

6.4- New directions 

Investigations of different masses, dimensionality, and quantum statistics are also 
producing results both for two (Barontini et al. 2009) and three-component systems 
(Williams et al. 2009, Wenz et al. 2009, Nakajima et al. 2011). Now more than one 
three-body system can be tuned to fulfill the Efimov conditions. Two components can 
form light-light-heavy and light-heavy-heavy combinations with very different scaling 
properties as described above. For three different particles the pair interactions involve 
three different scattering lengths and a correpondingly more complicated analysis. To 
be in the universal window at least two scattering lengths must be large compared to 
the effective ranges of the interactions. 

Extension to four-body universal states and corresponding recombination via four- 
body universal structures has been predicted (Hammer & Platter 2007, von Stecher 
et al. 2009) and now also experimentally observed for cold gasses (Pollack, Dries 
& Hulet 2009, Fcrlaino et al. 2009). The general structures of Brunnian systems 
(no bound subsystems) can be anticipated to be universal provided the states are 
located outside the effective range. Such investigations are only barely conceived at 
the moment in cold gases. In nuclei such systems exist but are too spatially confined to 
be truly universal (Curtis et al. 2008, Muta et al. 2011). One of the latest experiments 
and corresponding theoretical work even seem to suggest that five-body features could 
be accessible in atomic gases (Zenesini et al. 2012). 

A particularly surprising recent development has been the realization that the 
three-body parameter seems to be universal when expressed in units of the two- 
body interaction scale of the inter-atomic potential, which is the van der Waals 
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length, r y( ryy. The three-body parameter was denoted po in the coordinate-space 
hyperspherical approach discussed above but is also often quoted as a momentum- 
space cut-off (typically denoted A* or similar). Since po determines the lowest-lying 
Efimov state, the experimental findings can be written in terms of the threshold for 
the lowest Efimov state to appear out of the three-atom continuum on the a < 
side of a Feshbach resonance. We denote this threshold a_. The surprising result 
is that measurements on different atoms and using different Feshbach resonances 
yield a_/r yc j^ ~ —9.1 to within about 15% accuracy (Berninger et al. 2011, Wild 
et al. 2012, Knoop et al. 2012). This implies that there is some generic universality 
even in the three-body parameter which cannot be captured by simple zero-range 
models that need the po supplied from elsewhere. A number of theoretical works have 
presented various models that explain the observations (Naidon et al. 2012, Chin 2011, 
Wang, D'Incao, Esry & Greene 2012, Schmidt et al. 2012, Wang, Wang, D'Incao & 
Greene 2012) and it seems clear that the two-body inter-atomic potential is the culprit 
since it has a large repulsion at short distance which will naturally provide a three- 
body cut-off (S0rensen et al. 2012a, S0rensen et al. 20126). However, there is still a 
question of how the number of bound two-body states in the inter-atomic potential 
influences a_ (Wang, D'Incao, Esry & Greene 2012, S0rensen et al. 2012a). 

In the case of nuclear physics, it seems clear that the dense environment and the 
complicated nuclear interaction should not allow such an easy relation. For the sake 
of argument, imagine that we find an Efimov state in some nuclear system with a 
corresponding a_/r , with r the nuclear interaction range. If this number is of order 
-10 or so, then this implies that the nuclear system is universal since a_ is so much 
larger than r . However, as we have discussed at length in the previous sections, this 
is very unlikely except for the case of pure neutron matter. Unfortunately, as pointed 
out above, the three neutron system is unbound. 

Another interesting aspect in Efimov physics of recent times is the potential for 
Efimov states in the presence of long-range dipolar interactions. For fixed external 
polarization of the dipole moments (via applied magnetic or electric fields), it has 
been shown that the Efimov effect does occur for bosonic particles (Wang, D'Incao 
& Greene 2011a) while for fcrmionic atoms universal three-body states are possible 
(Wang, D'Incao & Greene 20116). This is interesting in comparison to nuclear 
physics (and condensed-matter physics) since long-range interactions like Coulomb 
are generically present. One could thus hope that the analogy can be used to 
learn something about nuclei. These dipolar system have also been studied in 
two-dimensional geometries (Wang et al. 2006, Wang 2007, Armstrong et al. 2010, 
Volosniev, Zinner, Fedorov, Jensen & Wunsch 2011, Zinner et al. 2012) and one- 
dimensional tubes (Klawunn et al. 2010, Wunsch et al. 2011, Zinner et al. 2011, Knap 
et al. 2012, Volosniev et al. 2013), where few-body bound states are even more prolific 
(Armstrong et al. 2012d, Volosniev et al. 20126). 

Efimov bound states also play a role in the recent exploration of the unitarity 
limit, where the so-called Tan relations can be used to deduce properties of 
bulk many-body systems from basic knowledge of few-body quantities (Tan 2008 a, 
Tan 20086, Tan 2008c, Braatcn & Platter 2008, Combescot et al. 2009, Barth & 
Zwerger 2011, Pricoupenko 2011, Valiente et al. 2011, Langmack et al. 2012, Valiente 
2012, Hofmann 2012, Valiente et al. 2012, Werner & Castin 2012a). It turns out 
that macroscopic observables such as the momentum distribution of a two-component 
Fermi system with \a\ — > oo are universal and depend on only one parameter called 
the contact, C (Stewart et al. 2010, Kuhnle et al. 2010). However, in the case 
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of bosonic systems that allow Efimov three-body states, there are in fact both a 
parameter for two- (C 2 ) and three-body contributions (C 3 ) (Braaten et al. 2011, Castin 
& Werner 2011, Werner &; Castin 20126). Even the fact that Efimov physics does not 
occur in two dimensions can be observed by measuring the contact parameters (Bellotti 
et al. 20136, Werner & Castin 20126). As the contact paramters are expected to be 
universal, they should be the same for a nuclear system in the limit of large scattering 
length. Measuring the momentum distribution during nuclear break-up could perhaps 
be a way to test predictions from this universal theory of strong interaction and 
compare to the measurements in cold atomic gases. This could in turn teach us 
something about the strongly-coupled nuclear environment. 

7. BCS-BEC in Nuclear Physics 

The basic physics of the BCS-BEC crossover is contained in the Fermi gas of a 
homogeneous system of particles as described in Section [4j The parameters are 
the Fermi momentum kp, or the density, and the s-wave scattering length a. In 
the ultracold gases one can fix the density n and vary the interaction parameter a 
through a Feshbach resonance. This variation takes the system from the BCS regime 
of unbound pairs, through the unitarity limit of | a | 1 = 0, and into the BEC regime 
of strongly bound pairs. The unitary regime can also be approached from either side 
by varying n for the Fermi gas since n|a| 3 (or fc^|a|) is the important parameter. This 
observation is the basis for attempts to find crossover effects in nuclear physics. Note, 
however, that it is not possible to cross unitarity by density variation and one will 
therefore be restricted to either BCS or BEC side of the crossover. 

7.1. From infinite matter to neutron halo nuclei 

As discussed earlier and listed in Tab. ([T]), the very large bare neutron-neutron 
scattering length implies that a large system of neutrons might reveal universal 
behavior corresponding to the limit 1/a = 0. Such an example is a neutron star, 
essentially infinite neutron matter, which perhaps might be described as the universal 
point in a BCS-BEC crossover model. 

Further investigations require an interaction but as we have discussed the bare 
nuclear force is only known phenomenologically, and the bare and effective in-medium 
nuclear interactions are different. This was evident in the discussion of polarization 
effects on the gap in Section [3] The details of the crossover strongly depends on the 
residual interaction beyond the mean-field approximation. Therefore to proceed we 
should separate mean-field and residual interactions. To illustrate, the BCS theory 
with the bare nuclear interaction produce a much smaller pairing gap in infinite 
nuclear matter than measured (Fetter & Walecka 1971, Dickhoff & Neck 2005). To 
improve one can introduce an effective and modified nuclear interaction with desired 
division between mean-field and residual interaction. However, this is then only 
applicable for a specific Hilbert space and a corresponding (BCS) approximation. 
The theoretical advantage of infinite nuclear matter is that the density and the 
corresponding interaction can be treated as a parameter. The crossover physics can 
then be studied as function of density. 

We now turn to realistic finite nuclear systems with wave functions obtained from 
appropriate nuclear interactions. We immediately face several fundamental issues 
concerned with validity or transfer of the concepts to such systems. The interaction is 
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in principle fixed and at best results for the universal point can be obtained. Choice of 
interaction and separation into mean-field and residual parts could perhaps allow some 
freedom but then the chosen input conditions directly control the output. The finite 
systems further introduces conceptual problems of coordinates, the spurious center- 
of-mass motion of a self-bound system, and the definition of condensates in terms of 
density matrices as discussed in Section ^ . We will focus on extracting the universal 
features related to the basic crossover concepts. Therefore we first need to address 
the seperation of mean-field and residual interaction. The troubles of nuclear physics 
in this respect are no less pronounced in nuclear studies of predominantly neutron 
matter. Here structures (such as a bound state of two neutrons, see below) can arise 
purely out of particular choices of residual pairing interaction between the nucleons. 
Suitably adjusted zero-range interactions to describe the residual interaction between 
neutrons in finite nuclei have been used for many years (Skyrme 1956, Bertsch & 
Esbensen 1991). On the surface of it such interactions appear similar to the zero- 
range or short-range approximations of atomic physics. There is, however, problems 
with this analogy as we now explain. 

The neutron-neutron interaction is often referred to as the pairing interaction 
which in nuclear physics is distinctly different from the full two-body interaction 
(Bertsch & Esbensen 1991). A commonly employed neutron-neutron pair residual 
interaction (Dobaczewski et al. 2001) has the form 



where r = |r x — r 2 | and R = |r x + r 2 |/2 are the relative and center-of-mass coordinates 
of the neutrons and p(R) the nuclear density at R with p the typical nuclear saturation 
density (Bohr & Mottelson 1975). The parameters ij and a determine the density- 
dependence of the interaction. Such density-dependent interactions are not used in 
atomic physics and any analogy must therefore be carefully examined. As one can 
see, for < r\ < 1 and < a < 1, the pairing strength in Eq. ( [49] ) will be peaked on 
the surface of the nucleus, something we will return to below. Unfortunately it is not 
clear what values are consistent with experiments (Dobaczewski et al. 2001), and the 
residual interaction is therefore somewhat arbitrary. 

Recent studies have made it clear that particularly the a parameter is crucial to 
get consistent, physically sound results (Dobaczewski et al. 2001, Rotival et al. 2009). 
The problem is that using r/ = 1 (strong pairing at the surface) and small a < 1 
can lead to crossover from BCS to BEC of neutron matter because the interaction 
unphysically binds the di-neutron system in the zero-density limit (Baldo et al. 1995). 
This leads, among other things, to prediction of neutron-rich halo nuclei that are much 
too large (Dobaczewski et al. 2001, Rotival et al. 2009). In particular, the pair density 
was found to have a constant profile out to large radii in Hartree-Fock-Bogoliubov 
calculations with a < 0.5 (Rotival et al. 2009) due to the bound di-neutron. A recent 
study using the BCS gap equations (Matsuo 2006) recommends potentials with a < 1 
for use in nuclear calculations, which is in conflict with the observations that this binds 
the di-neutron system. In the latter study, the features of the neutron wave function 
were seen as signatures of BCS-BEC crossover behavior. However, the troublesome 
potential means that this conclusion is unlikely to be transferable to other nuclear or 
atomic systems. 

The BCS-BEC crossover behavior seen in some nuclear studies using these 
parameters is then natural as the interaction favors it, but it might not be physical. 




(49) 
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However, such density-dependent interaction can still be used to reliably calculate 
observables such as energies (Bertsch & Esbensen 1991), and also the pairing gap in 
nuclear matter (Baldo et al. 1995). Thus whereas average quantities like energies are 



accessible, the density-dependent pairing interactions like Eq. (49 1 can be unphysical 



in their predictions for the spatial structure of the neutron pair wave function. 
7.2. Crossover in Finite Halo Nuclei 

We now discuss recent studies of BCS-BEC crossover physics in finite nuclei. A 
bound state of two neutrons in a potential from an ordinary core-nucleus was recently 
introduced as an example of BCS-BEC crossover (Hagino et al. 2007, Horiuchi & 
Suzuki 2006 a, Horiuchi & Suzuki 20066). Thus a bosonic state of two neutrons is 
suggested to display both BCS and BEC behavior. This is truly the simplest scenario 
imaginable. 

One of the most prominent examples of a nuclear three-body system is 11 Li 
( 9 Li+n+n) where the two neutrons are distributed far outside the core- radius R c . 
This is a Borromean halo system where none of the two-body subsystems are bound. 
Valence and core coordinates separate and the corresponding degrees of freedom 
decouple. The density variation then changes from normal nuclear density to an 
exponentially decreasing tail of the two neutrons. This interesting system was 
considered within BCS theory more than a decade ago (Barranco et al. 2001). 

The physics of the halo nucleus 11 Li has been discussed in the context of the BEC 
and BCS concepts. The analogy is deduced from the behavior of the three-body wave 
function, ^(r, R), described in one set of Jacobi coordinates, where r = r 1—1% connects 
the neutrons and R — (f\ + r 2 )/2 connect their center-of-mass to the core-nucleus. 
The s-wave part can then be projected out of \1/ and features of the remaining radial 
wave function fo(r, R) can be discussed in detail as function of the two coordinates 
(Hagino et al. 2007). The probability |/o| 2 exhibits an oscillatory structure as function 
of r for fixed small values of R below about 1.5 fm. These oscillations disappear when 
R becomes larger than about 4 fm where only one peak remains. This is qualitatively 
the BCS-BEC crossover behavior as function of the strength of an attractive residual 
interaction, i.e. first the BCS oscillations and then the BEC bound state. The other 
potential signature of crossover is the root mean square radius of fa with respect to r, 
again as function of R. This average distance between the two neutrons goes through 
a minimum roughly when R is equal to the core-radius. This is a qualitative feature 



of crossover since the extension of a pair- wave function, as in Eq. (26), is decreasing 
from the BCS to the BEC-side because the attraction becomes stronger and finally 
binds (Eagles 1969, Leggett 1980, Leggett 2006). 

We now discuss this interpretation of halo nuclei in terms of crossover physics 
to isolate the universality of the features of the physics involved and illustrate the 
influence of a particular model of the nuclear system under study and its interactions. 
The oscillatory behavior of the wave function and the minimum of the radius of the 
pair wave function must be generic, universal, and model- independent; unique signals 
only arising through the mechanism describing the BCS-BEC crossover. However, 
before discussing the question of universality, we caution that the suggested signatures 
involve non-observable quantities for the particular finite halo nuclei. It is therefore 
an indirect way to study crossover physics. 

The oscillating probability as function of r is a structure appearing through 
correlations of different single-particle orbits in BCS theory (Leggett 2006). To assess 
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Figure 13. The probabilities as functions of r for different values of R for 
two non-interacting particles in the mean-field from a square well potential with 
depth 108 MeV and radius 3 fm. Both particles are in the second s-state each 
with energy -0.17 MeV. Taken from (Zinner & Jensen 20086). 



whether this is a universal feature in a finite nucleus, we take W to be the uncorrelated 
product of two non-interacting single-particle neutron- 9 Li wave functions. We choose 



the second weakly bound s-wave in a square well potential and show in Fig. (13) 



the probabilities defined in (Hagino et al. 2007). The resemblance to the findings in 
(Hagino et al. 2007) is striking. However, the oscillations at small R are directly due 
to the nodes of the single-particle wave functions. They disappear when the lowest 
s-wave without nodes is used. As R increases these nodes are no longer geometrically 
compatible with at least one particle in the attractive part of the potential. Then 
the structure turns into one peak (Zinner & Jensen 20086). Adding a short-range 
neutron-neutron attraction of moderate strength would maintain these features and 
keep the neutrons a little closer resulting in a slower increase of the peak position with 
R for large R. 

We are led to conclude that the oscillatory behavior is not connected to many- 
body BCS-correlations, they are due to nodes of the single-particle neutron-core wave 
functions. Furthermore, the peak for large R is not a signal of a bound state or of 
condensation. This simple model reproduces the features which clearly demonstrates 
how careful one must be when assessing potential generic properties. It is always 
necessary to use more than one particular model of the system to check such behavior. 

The second feature, a minimum in the root mean square neutron-neutron distance 
at a finite center-of-mass distance, is interpreted as evidence for change of BCS- 
correlations into BEC-structure with R. The minimum itself can be related to the 
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choice of interactions in (Hagino et al. 2007), which closely resembles Eq. (49) with 
rj < 1 and a = 1. The neutron-core interaction is attractive inside the core radius R c 
and very small outside, and the neutron-neutron interaction is a J-function in r with 
an attractive strength which essentially is zero below R c and full strength for larger R. 
When R = the two neutrons can move independently inside the core, for R = R c the 
two neutrons can fully exploit their mutual attraction for r = 0, and for larger R they 
have to move apart to benefit from the more important neutron-core attraction. This 
produces the minimum in the rms value for R s» R c . We illustrate this in Fig. (14) 
by assuming a <5-shell potential at R = R c = 3 fin. The key properties are that the 
neutron-neutron attraction increases with R while the neutron-core attraction vanish 
outside the core. 
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Figure 14. The root mean square distance between the two neutrons as function 
of their center-of-mass distance from the core. Deep (dot-dashed) and shallow 
(dotted) gaussian neutron-core potentials with one and with no bound states are 
used. We also show the results when only s-waves are retained for the deep 
potential (dashed), and for non-interacting neutrons in a delta-shell potential 
around the core (full). 



Since n Li is very well studied in other successful three-body models, see e.g. 
(Garrido et al. 1999), we can compare results to test the universality requirement. 
We show in Fig. ( 14 ) neutron-neutron distances for a realistic bare neutron-neutron 
potential and different neutron-core potentials, i.e. a shallow and a deep gaussian 
without or with one bound state, respectively. These potentials are both able to 
reproduce the known three-body observables for Li. They differ at short distances 
where the neutron is inside the core but this inner part is unimportant in descriptions 
of spatially extended halo nuclei. The behavior seen in Fig. (14) differ significantly 
from that of (Hagino et al. 2007). Either no minimum when the lowest s-state is 
populated or an oscillation with one maximum and one minimum when the second 
s-state is populated. 

To see the effect of mixing partial waves we omitted p-waves, i.e. included only 
the s-waves in all Jacobi systems for the deep potential. The rms-curve remains similar 
but pushed to higher values. Finally we used a delta-shell neutron-core potential for 
two non-interacting neutrons. For a strongly attractive strength, where the neutrons 
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independently prefer to stay at R c , the rms-curve should go from 2R C through a small 
minimum at R = R c to a linearly increasing function at larger R. This only happens 
for very strong attraction otherwise the wave function is too smeared out although the 
same qualitative behavior remains. The conclusion is that generic pair correlations are 
not responsible for the rms minimum, which may appear for product wave functions 
and may be absent for correlated wave-mixed particles. 

We have to conclude that features of BCS-BEC crossover in two-neutron halo 
nuclei are not generic signals and the behavior is not universal. The example 
demonstrates the troubles that one faces in finite systems when comparing to BCS and 
BEC concepts. It also demonstrates some of the inherent dangers with specific residual 
nuclear interactions when they are pushed beyond their designated boundaries. 

8. Conclusions and perspective 

The status in the fields of cold atomic gases and nuclear physics differ at the moment 
very much from each other. The cold atomic gases have well controlled, simple 
interactions and obey simultaneously simple dilute limit conditions. Nuclei are dense 
self-bound many-body systems with complicated phenomenological interactions but 
well-controlled advanced techniques. Concepts and techniques from these fields can 
support each other. Theoretical methods developed for nuclei are very suitable for 
investigations of moderate numbers of bosons and fermions and mixtures of non- 
identical particles. Concepts originating from condensed matter and nuclear physics, 
now investigated in cold atoms may reappear in nuclei, and experimental techniques 
will soon allow investigations of systems with striking similarity to finite nuclei, but 
with controlled variation of the interactions. 

We have mainly focused on the simplest structures obtained in mean-field 
approximations. For bosons this implies that condensation could occur while for 
fermions pair correlations can lead to a superfluid state. In the fermionic case, 
the BCS-approximation is revealing after analogous degrees of freedom are identified. 
Pairs in time-reversed states in nuclei correspond to particles in two different internal 
(hypcrfine) states but the same set of external confinement quantum numbers (usually 
parabolic traps). Relative s-waves are by far the dominating contributors to the 
interaction in the zero-range approximation. We consider how to relate BCS and 
BEC coherent quantum states from condensed matter and cold atomic gas physics to 
states in finite systems of few particles. The goal is to discuss universal features of 
these quantum states and phases and investigate whether they can be manifest in few- 
body systems with nuclei as the example. A particularly interesting venture at a time 
when trapping of small samples of only a few atoms have been realized experimentally 
(Serwane et al. 2011, Ziirn et al. 2012, Bakr et al. 2010, Shcrson et al. 2010, Weitenberg 
et al. 2011), implying that systems resembling few-body nuclei can be produced and 
studied in cold atomic gase experiments. 

An important contrast is the fact that cold atomic gases are held together by 
an external field. For repulsive two-body interactions, an almost perfect BEC is 
possible (Anderson et al. 1995, Davis et al. 1995). The influence of correlations 
are small because the system is extremely dilute (Pethick & Smith 2008) (away 
from two-body scattering resonances, i.e. Feshbach resonances). For moderate 
attractive interactions and relatively few particles similar stable structures also exist 
(Ruprecht et al. 1995, Baym & Pethick 1996). However, attraction destabilizes these 
systems and collapse occurs into very dense and highly correlated structures (Bradley 



CONTENTS 



60 



et al. 1997, Roberts ct al. 2001, Donley et al. 2001, Lahaye et al. 2008). Nuclear 
systems become similar to large samples of cold atoms only in the dilute limit when 
the density is about two orders less than that of nuclear saturation. Attempts to 
approach the limit of dilute nuclear states are found in the study of cluster states or 
weakly bound states close to particle or cluster thresholds. The BEC analog has been 
introduced as a few, or many, a-particles in a dilute state. However, the Coulomb 
repulsion would lead to immediate fragmentation or, if binding is achieved through 
the attractive nuclear interaction, the particles must be close-lying as for ordinary 
nuclei. This in turn implies that the diluteness criteria is violated and the BEC 
feature destroyed. Furthermore, detecting a condensate requires knowledge of the 
coherence in the system, for example through a two-particle correlation function. In 
the case of nuclear a-condensates this is impossible to measure and only standard one- 
body density information can be obtained. The notion of condensate in finite nuclei 
is therefore ill-defined and should be regarded as a crude approximation to the cluster 
structure of excited states in a-nuclei. Attempts to avoid the Coulomb interaction by 
using neutrons suffer from a similar deficiency, i.e. two neutrons can be weakly bound 
to a core and form a neutron halo, but adding more than a few extra nucleons would 
either make the system unstable due to the Pauli principle or lead to ordinary bound 
nuclear states. Many-body neutron halo states would fall apart or form clusters of 
two, three or more pieces. The external confining field is not present. Combinations 
allowing bosons as deuterons does not seem promising for the same reasons, i.e. the 
system becomes either unstable due to lack of attraction or the spatially extended 
dcuteron structure is lost as soon as nuclear surface densities arc reached. 

A promising direction of cross-fertilization between nuclear physics and cold 
atomic gases has turned out to be the investigation of three-body bound states in the 
universal low-energy regime. The Efimov effect where an infinite number of three-body 
bound states appear near the threshold for two-body binding was initially predicted in 
the context of nuclear physics. With its recent experimental observation in cold atomic 
gases it has seen a great revival and tremendous activity in the area has ensued. Here 
we have discussed the universal predictions from the original hyperspherical point of 
view, and have suggested systems which have favourable parameters for observing the 
effect. The flexible tunability of experimental parameters in cold atomic gases can 
help elucidate the features of the Efimov effect, also beyond the universal predictions. 
This information could then be transfered to nuclear physics and aid in understanding 
why it is so difficult to see the Efimov effect in nuclear systems. 

Understanding cold atomic Fermi gases are to a large extend based on pairing 
Hamiltonians, and the related BCS-type approximations. This is a much better 
approximation than for nuclei, but the theoretical models often employed in atomic 
physics are not always as accurate as the impressive experimental precision. In 
nuclear physics, improved models had to be developed, and these have now proven 
their capability in the description of many-fermion systems. Examples are large-scale 
shell-model methods, hyperspherical basis expansions, algebraic methods for pairing 
Hamiltonians, and density functional theory. They are based on first principle and go 
far beyond mean-field and pairing correlations. They are highly flexible in connection 
with choice of interactions and should thus be able to provide accurate results also 
for atomic gases. The methods are based on effective interactions adapted to reduced 
Hilbert spaces which is a tremendous advantage for the atomic A-body systems. These 
methods were first of all developed for handling identical fermions (nucleons), but 
also simpler bosonic systems can in most cases be treated with the same numerical 



CONTENTS 



61 



techniques. Strongly-coupled systems at unitarity and the BCS-BEC crossover region 
are thus obvious fields to explore. Also finite-range effects and other more complicated 
corrections are naturally handled by nuclear methods. 

In summary, we have discussed the physics of BEC, BCS, and BCS-BEC crossover 
as studied in cold atomic gases from a nuclear physics perspective. We emphasized 
similarities, and differences through a detailed discussion of some relevant methods in 
wide use. We find that finite size effects in nuclear systems will most often destroy 
signatures of the coherent behavior that is routinely observed in cold atomic gases. No 
current models or measurements seem to support the interpretation of nuclear states 
as manifestations of BEC, BCS or crossover physics (with one potential exception 
being infinite neutron matter). At the moment it appears to be difficult to specify 
distinguishing measurable consequences or to design test experiments of the required 
generic features of universality, a basic first step in any successful transfer of these 
concepts to nuclei in particular and few-body systems in general. However, as cold 
atomic gas experiments with few particles become reality the situation is likely to 
change. On the other hand, here we have suggested some tempting cold atom systems 
to investigate with nuclear techniques in order to facilitate further exchange between 
the fields. 
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